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Summary 


This  paper  contains  an  asymptotic  investigation  -  within  the  nonlinear 
theory  of  elastostatic  plane  stress  -  of  the  deformations  and  stresses  near 
the  tips  of  a  traction-free  interface-crack  between  two  dissimilar  semi-infinite 
Neo-Hookean  sheets.  The  results  obtained  are  free  of  oscillatory  singularities 
of  the  kind  predicted  by  the  linearized  theory,  which  would  require  the  two 
deformed  faces  of  an  interface-crack  to  overlap  in  the  vicinity  of  its  tips. 
Instead,  the  crack  is  found  to  open  smoothly  near  its  ends,  regardless  of  the 
specific  loading  at  infinity. 


1 .  Introduction 

Owing  to  its  importance  in  fracture  mechanics,  the  plane  problem  of  one 
or  more  interface-cracks  between  two  dissimilar  elastic  slabs  has  received 
repeated  attention  in  linearized  elastostatics.1  The  earliest  such  investiga¬ 
tion  appears  to  be  due  to  Williams  [1]  (1959),  who  examined  the  local  character 
near  a  tip  of  a  traction-free  interface-crack  of  elastostatic  fields  compatible 


* 

The  results  communicated  in  this  paper  were  obtained  in  the  course  of  an  in¬ 
vestigation  supported  in  part  by  Contract  NOOOT 4-75-C-01 96  with  the  Office  of 
Naval  Research  in  Washington,  D.C. 

complete  bibl " •'qraphy  of  the  extensive  literature  on  this  subject  is  beyond 
the  scope  of  thi  aper.  Additional  references  can  be  found  in  those  listed 
at  the  end  of  the  paper. 
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wlth  the  governing  field  equations,  as  well  as  with  the  appropriate  boundary 
and  bond  conditions.  The  analysis  in  [1]  is  based  on  an  asymptotic  scheme, 
originated  by  Kneln  [2],  which  Williams  [3]  had  employed  previously  to  explore 
two-dimensional  corner  singularities  possible  in  homogeneous  elastic  slabs. 

The  scheme  here  alluded  to  leads  from  the  biharmonic  eauation  satisfied  by  the 
generating  Airy  stress  function  to  an  eigenvalue  problem  for  an  ordinary  fourth- 
order  differential  equation.  In  [1]  the  prevailing  boundary  and  continuity 
conditions  give  rise  to  a  sequence  of  complex  eigenvalues  that  induces  oscil¬ 
latory  singularities  at  a  tip  of  an  interface-crack  in  the  associated  sequence 
of  displacement  and  stress  fields.1 

A  global  solution  to  the  problem  of  a  uniformly  pressurized  interface- 
crack  of  finite  length,  between  two  homogeneous  and  isotropic  (linearly  elastic) 
semi-infinite  slabs  of  possibly  distinct  material  properties,  was  deduced  by 
England  [4]  (1965).  Let  2z  be  the  length  of  the  crack  and  choose  rectangular 
Cartesian  coordinates  (x-pX2)  as  indicated  in  Fia.l.  Next,  let  Hj  and  H2 
be  the  half-planes  x2  >  0  and  x2<0,  occupied  by  the  interior  of  the  middle 
cross-section  of  the  upper  and  lower  slab,  respectively,  while  denoting  by 
and  ( k  =  1,2)  the  corresponding  shear  moduli  and  Poisson-ratios  of  the  two 

O 

slabs.  If  u  and  a  „  are  the  components  of  displacement  and  stress  in  the 
a  aB 

underlying  coordinate  frame,  England's  solution  conforms  to  the  followino  formu¬ 
lation  of  his  problem:  it  satisfies  the  two-dimensional  displacement-stress 
relations  and  stress  equations  of  equilibrium  (in  the  absence  of  body  forces), 

1  Similar  oscillatory  singularities  -  likewise  traceable  to  complex  eigenvalues  - 
are  encountered  in  [3]  in  connection  with  certain  mixed  boundary  conditions. 

2From  here  on  the  "material  index"  k,  as  well  as  all  Greek  subscripts,  are 
understood  to  take  on  the  values  (1,2);  further,  Greek  subscripts  are  subject 
to  the  usual  summation  and  differentiation  conventions. 
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so  that 


‘afT^k 


Kk  "  3 

1  -  5a8Uo  ,(j  +  l  U(a,e) 


•V.fT0  on  Hk  • 


(l.l)1 


where  6  „  is  the  Kronecker  delta  and 
ae 


*  3  -  4\>k  for  plane  strain, 


3  “  v  k 

Kk  =  T+~7"  ^or  9enera^ize£i  plane  stress  ; 


(1.2) 


it  obeys  the  loading  conditions 

cr2e^xl’*0±^  =  -<S2sa*  (0>°)’  C I  x1 1  <  0 

2  2 

oae(xi  ,x2)  =  °(1)  as  x1+.^co  , 


-\ 


> 


j 


(1.3) 


in  which  a  is  the  constant  pressure  applied  to  either  crack-face,  and  it 
meets  the  bond  conditions 


ua(x1 ,0+)  =  u(j(x1 ,0-),  a2e(x1 ,0+)  =  o^x,  ,0-)  (|x1|>t)  .  (1.4) 


In  addition,  u  is  twice  continuously  differentiable  on  H.  ,  the  limits 

a  K 

u  (x, ,0±)  exist  and  depend  continuously  on  x,  for  -oo  <  x,  < oo  ,  whereas 

01  I  II 

the  limits  c  „(x,,0±)  exist  and  are  continuous  for  all  x,  f  ±£. 

dp  f  I 

It  should  be  recalled  that  the  generalized  plane-stress  solution  pertains 
to  "thin"  slabs  (elastic  sheets).  In  this  instance  u  (x,  ,>u),  a  ,(x,,x_) 

a  t  c  ap  I  c 

at  once  approximate  the  thickness-averages  and  the  mid-plane  values  of  the 
corresponding  displacements  and  stresses. 

^Here  uj  7,  =  (u  7+u„  )/2  are  the  infinitesimal  strains. 

(a, 6)  a, 8  6, a 
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We  now  introduce  local  polar  coordinates  (r,e)  with  the  oriqin  at 
x2  *  0  (see  Fig,  1 )  and  describe  the  asympotic  behavior  near  the  right- 
hand  crack-tip  of  the  solution  obtained  in  [4],  If  this  tip  is  fixed  by  super¬ 
position  of  a  suitable  rigid  translation  upon  the  displacement  field  at  hand, 
one  has  as  r  +  0, 


ua~/r  [cfa(e)cos(Y  log  r)  +  ijJa(e)sin(y  log  r)] 


oa6~  -jr[cfa6(e)cos(Ylogr)+^ae(8)sin(Ylogr)]  , 


0.5) 


where  cp^,  ^  are  fully  determinate  functions  of  the  polar  angle 

9,  and  v  is  the  material  parameter  defined  by 


1  s  +  K.j 

y *  27  l09a*  a*  s*  7J  • 


(1.6) 


s  being  the  shear  stiffness-ratio  of  the  two  slabs.  Equations  (1.5)  bring 
into  evidence  the  unboundedness  of  o  at  the  crack-tips,  as  well  as  the 

0t(5 

oscillatory  nature  of  the  resulting  crack-tip  singularities. 

With  a  view  toward  examining  the  difficulty  that  led  England  to  reject 
his  solution  as  physically  inadmissible,  we  cite  from  [4]  the  formula  for  the 
relative  normal  displacement  of  the  two  crack-faces: 


a2(x] )  =  u?(x1 ,0+)  -  u2(x1 ,0-)  = 


a/a 

2(1  +a) 


r~  1  +  <i  1  Ko  1  /  o  9  /  i  +  X,  \ 

*xi  C0S(T,09r^r;-)  (|xi'<il)  •  <1-7> 


When  y^O*  (1.7)  reveals  an  infinity  of  oscillations  in  the  sign  of  a2(x^) 
as  X-J  +  +&  and  thus  implies  wrinkling  of  the  faces  of  the  crack  in  the 
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vicinity  of  its  ends.^  In  order  to  infer  the  unacceptable  prediction  of 
over  1  appi ng  crack-faces  it  is  essential  to  consider  also  the  corresponding 
relative  tangential  displacement,  which  is  readily  computed  with  the  aid  of 
results  given  in  [4]: 


(xi )  -  u1(x1 ,0+)  -  u1(x1 ,0-)  = 


2TH^T 


^K1  ,  U<2 


J1 


/t2-x^  sinolog  (|x]!<0  .  (1 


8) 


Equations  (1.7),  (1.8)  evidently  imply  the  presence  of  interpenetration  if 
there  is  a  value  of  satisfying 

^  (x-| )  =  0,  a2(x-| )  <  0,  jxj  |  <  s.  •  (1.9) 


If  y  /  0,  one  confirms  easily  that  (1.9)  in  fact  possess  a  doubly  infinite 
sequence  q*  roots  supplied  by 


xj  j  ^  =  4  tanh( jff/2y)  (j  =  ±1  ,±3,±5, . . . ) 


(1.10) 


From  (1.10)  follows 


<5  =  max  [4  -  |xjJ )  |j  =  4[1  -  tanh(n/2iy| )]  (1.11) 


•8 

and  one  finds  that  6/4  <4x10  for  all  physically  realistic  values  of  the 
2 

elastic  constants.  Consequently,  the  roots  of  (1.9)  are  confined  to 


Since  uk  >  0,  -l<vjcil/2,  it  is  clear  from  (1.2),  (1.6)  that  the  material 
parameters  ,  s,  and  a  are  positive;  li<k<7  for  plane  strain,  while 
5/3*<|<<oo  for  generalized  plane  stress.  Thus  A?(xi)  cannot  vanish  iden¬ 
tically. 

2uk  >0,  Oiv^il/2. 
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exceedingly  small  intervals  adjacent  to  the  ends  of  the  crack.  Nevertheless 
the  presence  of  arv^  overlap  violates  the  physical  requirement  that  the  re¬ 
sulting  deformation  mapping  be  one-to-one  and  thus  -  as  pointed  out  in  [4]  - 
Invalidates  the  solution. 

Evidently,  the  foregoing  violation  disappears  if  and  only  If  y-0  or, 
equivalently. 


s  +  <1  =  1  +  sko  (s  =  u-,/u9)  , 


(1.12) 


In  which  instance  the  ensu'ng  singularities  are  no  longer  oscillatory  and 


A.j(x^)  —  0  ,  ^(x^ )  >  0  ( I  |  *  & )  . 


(1.13) 


Condition  (1.12)  is  met  in  particular  for  plane  strain,  regardless  of  the 
value  of  s.  In  the  limiting  case  of  two  incompressible  slabs,  since 
\»!  =  \>2  =  l/2  here  Implies  *-|  =  *  In  contrast»  1*  is  essential  for  our 

purposes  to  observe  that  (1.12)  fails  to  hold  for  generalized  plane  stress 
of  two  incompressible  sheets,  unless  the  two  materials  are  identical.^ 

We  turn  next  to  the  analogous  problem  for  a  traction-free  interface- 
crack  between  two  semi-infinite  slabs  that  are  subjected  to  loads  at  infinity. 
In  this  connection  we  confine  our  attention  to  the  case  in  which  the  loading 
conditions  (1.3)  give  way  to 


°2g(*i*®-)*®  ( (x-j  |  <  t )  , 

0n(xrx2)’SoiV  +  o(l)  as  x-j  +  Xj  -*■  oo,  (x-j  ^2)  £ 

°2b^X1  ,X2^  *  S2b°22  +  ®  ^  x.j  +  Xg -*■  00  ^22  ^ 


(1.14) 


Note  from  the  second  of  (1.2)  that  v|(=l/2  now  gives  <^=5/3,  so  that  (1.12) 
demands  y-,  ~  ^2' 
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As  was  made  clear  by  Rice  and  Sih  [5],  an  arbitrary  assignment  of  the  constants 

go  (  K  1  00  ,  00  . 

o^i'  and  egg  is  in  general  inconsistent  with  the  bond  conditions 

(1.4),  the  first  of  which  demands  the  continuity  across  the  bonded  interface  of 
the  extensional  strain  u-j  ^ .  In  view  of  the  displacement-stress  relations 
(1.1)  this  requirement  leads  to  the  loading  constraint 

(1  +  ic.j  )ojj  ^  -  s ( 1  +  ^2 )°1  ? ^  +  [*1  "  3  -  s(i<2  -  3)]?22  -  0  •  (1.15)^ 

Observe  that  o|^  =  0  is  admitted  by  (1.15)  if  and  only  if 

«1  -  3  =  s ( i-g  -  3),  ($  =  1^/^)  •  (1.16) 


Accordingly,  the  problem  of  an  interface-crack  in  a  uni-axial  tension  field 
(at  right  angles  to  the  interface)  cannot  possibly  possess  a  solution  unless 
(1.16)  holds.2 

When  the  given  load  parameters  o|^»  a22  obey  (1.15),  the  solution  to 
the  problem  governed  by  (1.1),  (1.4),  (1.14)  is  obtained  at  or.ce  by  superposition 
upon  England's  [4]  solution  for  the  uniformly  pressurized  crack  with  0  =  ^2 
of  the  piecewise  homogeneous  deformation  with  the  displacements  and  stresses: 


U1  =  [*Kk  +  1Ml^  +  (<k"  3)o22xl]  ’ 

U2S  5p7  [^lck  +  1^°22+  *Kk"  3^°1 1  ^2]  ’ 

k 


°n 


“(k) 

°11  ’  °22 


°22’ 


°12  _  °21  =  C 


on 


>  (1.17) 


Hhis  relation  is  equivalent  to  Eq.(19)  in  [5].  No  additional  constraints  on 
the  load  parameters  arise  if  ooi  "tends  to  a  constant  non-zero  value 

o12  at  infinity. 

p 

This  fact  appears  to  go  unnoticed  in  some  papers.  Note  that  (1.16)  is  satisfied 
not  only  In  the  special  case  s = 1 ,  v,  = (identical  materials),  but  also  for 
vk*0  (tck=3),  regardless  of  the  value  of  s. 
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Slnce  u  in  (1.17)  is  continuous  across  the  entire  x,-axis,  it  is 

a  -  I 

clear  that  (1.7),  (1.8)  and  the  discussion  followina  these  equations  remain 
valid  also  for  the  problem  of  the  traction-free  crack  under  present  consid¬ 
eration.  Further,  the  asymptotic  results  (1.5)  evidently  continue  to  hold 
In  the  present  circumstances.  Indeed,  the  right-hand  members  In  (1.5)  are 
consistent  with  the  displacements  and  stresses  appropriate  to  the  only  member 
of  the  sequence  of  elastostatic  fields  deduced  by  Williams  [1]  that  has  con¬ 
tinuous  displacements  but  unbounded  stresses  at  the  crack-tip.1  This  field 
encompasses  at  the  same  time  the  asymptotic  structure  near  a  crack-tip  of  the 
global  solution  associated  with  an  all-around  uniform  shear  loading  at  in- 

2  3 

finity,  which  is  once  again  described  by  (1.5). 

4 

The  objections  raised  in  [4]  to  violations  of  the  impenetrability  re¬ 
quirement  apply  to  a  host  of  papers  on  various  interface-crack  problems  that 
have  appeared  in  the  literature.  These  misgivings  have  prompted  a  renewed 
theoretical  concern  with  Interface-cracks  In  recent  years.  Thus  Comninou 
[7]^  (1977),  [8]  (1978)  sought  to  remove  the  Inadmissible  field  oscillations 
on  the  basis  of  the  ad-hoc  assumption  th>t  the  crack-faces  remain  in  friction¬ 
less  contact  over  two  sub-segments  (adjacent  to  the  two  ends  of  the  crack) 
of  initially  undetermined  length.  Still  more  recently  Achenbach,  Keer,  Khetan, 
and  Chen  [9]  (1979)  relied  on  the  Dugdale-Barenblatt  [10],  [11]  model  of 

'it  should  be  recalled  that  the  results  In  [1]  involve  certain  undetermined 
amplitude  coefficients  which  elude  tht  local  analysis  carried  out  there. 

2 

See  [5],  where  such  a  loading  is  included. 

3In  this  Instance  cpa,  <j<a,  cpag,  ipag  are  no  longer  the  same  functions  of  the 
polar  angle  as  before. 

4The  same  criticism  was  voiced  independently  by  Malyshev  and  Salganik  [6], 

^This  publication  contains  a  fairly  extensive  list  of  references  to  earlier  work. 


inelastic  behavior  near  the  tips  of  a  crack  in  order  to  eliminate  unbounded 
stress  singularities  altogether  from  problems  involving  interface-cracks. 


In  the  current  paper  we  adhere  to  strictly  elastic  behavior  but  relinquish 
the  hypothesis  of  infinitesimal  deformations,  which  is  in  fact  violated  by  solu¬ 
tions  exhibiting  locally  unbounded  stresses  and  hence  unbounded  displacement 
gradients.  We  aim  to  show  that  the  offensive  oscillatory  sinaularities  arisinn 
in  interface-crack  problems  stem  from  the  linearization  of  such  problems,  as 
conjectured  by  Ergland  [4],  rather  than  from  the  assumption  of  perfect  elasticity 
or  the  particular  idealization  underlying  the  formulation  of  the  boundary  and 
bond  conditions.  With  a  view  toward  accomplishing  this  purpose  in  an  analyti¬ 
cally  amenable  setting  we  deal  asymptotically  with  the  elastostatic  field  near 
a  tip  of  a  traction-free  crack  between  two  otherwise  bonded  incompressible 
slabs  of  Neo-Hookean  material.  Furthermore,  the  present  study  is  carried  out 
within  the  nonlinear  theory  of  elastostatic  plane  stress  since  a  plane-strain 
analysis  appropriate  to  incompressible  slabs  would  be  pointless.*  To  be  sure, 
the  theory  of  plane  stress  -  in  contrast  to  plane  strain  -  is  approximate  and 
presupposes  the  slab-thickness  to  be  small  compared  to  a  characteristic  in-plane 
dimension  of  its  cross-section.  Objections  based  upon  the  approximate  nature 
of  this  theory,  which  might  be  raised  especially  in  connection  with  crack  prob¬ 
lems,  however  apply  equally  to  linearized  and  finite  elastostatics  and  do  not 
interfere  with  our  primary  purpose. 

The  method  used  in  pursuing  the  local  issue  at  hand  is  an  adaptation  to 
the  nonlinear  theory  of  the  asymptotic  scheme  underlying  [1],  [2],  [3].  The 
same  approach  was  employed  in  a  sequence  of  previous  studies  of  crack-problems 

^Recall  from  the  discussion  following  (1.13)  that  the  1  inear  theory  does  not 
lead  to  oscillatory  singularities  in  these  circumstances. 
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in  the  nonlinear  equilibrium  theory,  under  various  loading  conditions  and 
diverse  constitutive  assumptions.  These  papers  are  referenced  and  briefly 
summarized  in  two  recent  survey  articles  [12],  [13],  Particularly  pertinent 
to  the  present  study  is  a  ^elated  asymptotic  exploration  [14]  of  elastostatic 
singularities  induced  by  certain  mixed  boundary  conditions  of  the  kind  arising 
in  the  problem  of  the  "rough  punch".  While  the  linearized  theory  in  these 
circumstances  predicts  oscillatory  singularities,  their  absence  is  shown  in 

[14]  to  be  consistent  with  the  nonlinear  theory  of  plane  strain  for  compres¬ 
sible  materials  of  the  harmonic  type. 

The  present  work  is  also  closely  related  to  a  paper  by  Wong  and  Shield 

[15]  (1969)  that  predates  the  publications  mentioned  above.  In  [15]  an  ap¬ 
proximate  global  plane-stress  solution  is  deduced  for  the  problem  occasioned 
by  a  finite  crack  in  an  all-around  infinite  Neo-Hookean  sheet,  subjected  to 
bi-axial  tension  at  infinity.  The  approximative  approach  adopted  there  re¬ 
quires  the  deformations  to  be  large  throughout  the  sheet.  Our  local  results 
pertaining  to  an  interface-crack  between  two  distinct  Neo-Hookean  sheets, 
upon  proper  specialization,  are  found  to  be  in  asymptotic  agreement  with  the 
solution  reported  in  [15]. 

The  chief  conclusion  reached  by  us  concerns  the  absence  of  oscillatory 
singularities  in  the  present  nonlinear  setting  of  the  interface-crack  problem. 
Somewhat  surprisingly  we  find  that  a  result  arrived  at  in  [15],  according 
to  which  the  crack  opens  up  smoothly  under  the  symmetric  loading  considered 
there,  continues  to  hold  true  for  the  interface-crack  -  regardless  of  the 
particular  loading  conditions.1 
— 

We  exclude  degenerate  loadings  that  give  rise  to  finite  displacement  gradients 
at  the  crack-tips,  such  as  uni -axial  tension  parallel  to  the  crack-faces  and 
compatible  with  (1.15). 
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Finally,  as  a  by-product,  the  analysis  presented  in  this  paper  yields 
certain  results  of  interest  in  connection  with  the  Mode  II  crack  problem  for 
a  single  (homogeneous)  Neo-Hookean  sheet.  The  global  solution  of  this  problem 
within  the  infinitesimal  theory  predicts  that  the  crack-faces  fail  to  separate 
in  this  instance  and  that  a  Mode  II  loading  of  simple  shear  at  infinity  leads 
to  an  elastostatic  field  which  is  anti-syirmetric  with  respect  to  the  crack- 
axis.  Both  of  these  predictions  are  found  to  be  in  conflict  with  the  local 
results  established  in  this  paper  and  hence  reflect  degeneracies  stemming  from 
the  linearization  of  the  problem.  Analogous  conclusions  regarding  the  Mode  II 
crack  problem  were  arrived  at  previously  by  Stephenson  [16]  in  a  nonlinear 
asymptotic  study  encompassing  a  class  of  incompressible  materials  under  con¬ 
ditions  of  plane  strain J  In  particular,  Stephenson  succeeded  in  provinci  that 
the  global  nonlinear  Mode  II  crack  problem  for  Diane  strain  of  a  Neo-Hookean 
material  cannot  admit  a  solution  anti -symmetric  about  the  crack-axis. 

2.  Preliminaries  from  the  theory  of  plane  stress  in  finite  elastostatics. 

Neo-Hookean  materials. 

As  prerequisites  for  the  analysis  to  follow,  we  assemble  in  the  present, 
largely  expository,  section,  some  basic  ingredients  of  the  nonlinear  equllib- 
rium  theory  of  plane  stress  and  in  this  connection  confine  our  attention  to 
incompressible,  homogeneous  and  isotropic,  elastic  solids. 

Let  (x.j,X2,x3)  be  rectangular  Cartesian  material  coordinates  and  con¬ 
sider  a  body  which,  in  an  undeformed  configuration,  occupies  the  closed 

See  also  [17],  where  nonlinear  effects  bearing  on  the  question  as  to  whether 
or  not  a  crack  opens  up  In  the  presence  of  a  Mode  II  loading  are  explored  further. 

2 

An  account  of  this  approximate  theory,  couched  in  general  tensor  notation,  is 
contained  In  a  paper  by  Adkins,  Green,  and  Nicholas  [18]  (1954),  who  cite  per¬ 
tinent  earlier  work.  See  also  [15]. 


cylindrical  region  ft  described  by 


ft=  (xj  (x^  ,x2)  €  n,  -t/2  s.  x3  s  t/2]  , 


(2. 1)1 


where  ji  is  the  cross-section  of  ft  in  the  mid-plane  -'<3  =  0  and  t  the 
constant  thickness  of  the  cylindrical  slab. 

A  locally  volume  preserving  deformation  of  the  body  is  characterized  by 

y(x)  =  x  +  u(x)  for  all  x6ft  ,  (2.2)^ 


with 


det  F  =  1 ,  F  =  vy  on  ft  , 


(2.3) 


in  which  and  F  stand  for  the  displacement  vector  field  and  the  deforma¬ 
tion-gradient  tensor  field,  respectively.  We  shall  temporarily  take  for 
granted  that  the  mapping  y  is  twice  continuously  differentiable  and  uniquely 
invertible  on  ft.  If  a  Is  the  nominal  (Piola)  stress  field  accompanying 
the  deformation,  equilibrium  -  in  the  absence  of  body  forces  —  demands  that 


div 


on  ft 


(2.4) 


Further,  a  Is  linked  to  the  corresponding  actual  (Cauchy)  stress  field  t 
by  means  of 


tsoF^  on  ft 

r*>  ^v/v 


( 2 . 5 ) 3 


For  the  time  being,  letters  in  boldface  denote  vectors  and  second-order  tensors 
in  three  dimensions. 

2 

Thus  yj(x,)  (1  *1,2,3)  are  the  spatial  coordinates  of  the  material  point 
x =  (x^  .Xg.x^)  after  the  deformation. 

3A  superscript  T  Indicates  transposition.  Note  that  at  present  ^  is  re¬ 
garded  as  a  function  of  position  on  R.  " 
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Next,  call  G  the  left  Cauchy-Green  deformation  tensor  associated  with  the 
deformation  (2.2)  and  let  I- ,  I3  designate  the  fundamental  scalar  in¬ 
variants  of  G.  Thus, 


G  =  F  F  ort  ft  , 


I-j  strG,  r2*  Jc(tr6)2- trG2],  I3  =  detG  =  l 


A 


r  (2.6) 


J 


Suppose  now  the  body  at  hand  possesses  an  elastic  potential  W(I^,I-)t  so 
that  W^,,^)  represents  the  strain-energy  density  per  unit  undeformed 
volume.  The  appropriate  constitutive  law  then  takes  the  equivalent  alterna¬ 
tive  forms 


r  =  AG-BG2-p1,  ct  =  AF- BSF- pF‘T  , 

/v  rw  /w  «v  r/v 


(2.7) 


where 

A  =  2(Wt  +  I.W.  ),  B  =  2W.  ,  (2.8)2 

X1  1  l2  lZ 

while  p  stands  for  the  arbitrary  scalar  pressure  needed  to  accommodate  the 
kinematical  constraint  of  incompressibility. 

Finally,  we  recall  that  if  x,.  (i  =  1,2,3)  are  the  local  principal 
stretches  associated  with  the  deformation  (2.2),  their  souares  are  the  local 
eigenvalues  of  the  symmetric  positive-definite  tensor  field  G.  Accordingly, 
(2.6)  gives 

I^x^  +  x^  +  x2,  I2  =  x2x2  + x2x2  + x2x2,  I3-*i*i>3“l  •  (2.9) 

^ere  1  is  the  idem  tensor  and  F’"'"  the  transposed  inverse  of  F. 

/v  /v  rw 

2 

We  write  and  for  the  partial  derivatives  of  W  with  respect 

to  the  corresponding  deformation  invariants. 


At  this  stage  we  subject  the  deformation  (2.2)  to  the  following  special¬ 
izing  assumptions: 

(a)  it  is  symmetric  about  the  mid-plane  x3  =  0,  so  that  ya(x-|  »x2>x3) 

is  an  even,  and  an  odd  function,  of  x^; 

(b)  it  is  consistent  with  the  requirement  that  the  deformed  faces  of 
the  cylindrical  slab  be  free  of  tractions,  whence 

ai3(xl  ,x2,±t/2)  =  0,  (xrx2)en,  (i  =  1,2,3)  .  (2.10) 

In  addition,  we  suppose  the  cylindrical  slab  to  be  "thin”  in  the  sense  that 
t  is  small  compared  to  a  characteristic  dimension  of  its  mid-section  n. 
Motivated  by  this  supposition,  as  well  as  by  (a)  and  by  (2.10)  with  i  =  3, 
we  adjoin  to  (a)  and  (b)  the  approximative  assumption 

(c)  °33'®  on  n’  °ote  3  =  ®  00  R  '  (2-11) 

O 

If  f  is  a  function  defined  on  R,  we  shall  write  f  for  the  restric¬ 
tion  of  f  to  n.  The  relevant  theory  of  plane  stress  alms  at  an  approxi¬ 
mate  characterization  of  y  ,  a  .,  and  °  (in  the  present  constitutive 

ci  dp  dp 

setting)  on  the  basis  of  assumptions  (a),  (b),  (c). 

In  view  of  (2.2),  one  has 

y  (x,,x9)  =  x  +u  (x,,x,)  on  n  .  (2.12) 

-'a  1  C  a  a  I  c. 

Further,  the  equilibrium  equation  (2.4)  implies 
t/2 

Ccla6,e(xl,x2,x3)+Oa3,3(xl’x2,x3^dx3=0  ’  (2J3) 

-t/2 

which,  because  of  (2.10)  and  the  second  of  (2.11),  reduces  to 
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1 


II 


°«s,6’0  on  "  ■ 


(2.14) 


Next,  upon  putting 


x(x1 ,x2)  =  y3j3(x1 ,x2,0)  on  n 


(2.15) 


(2.16) 


one  draws  from  (2.3),  (2.5),  (2.6),  (2.7),  together  with  (a),  that 


f«3-f3.-°-  ?33**'  i343„  =  °-  a33sl2  • 


°3a  =  So3  *  0  ’  “33  =  T(*x2  '  ‘  ”p)  • 


Ta3  "  T3a  ~  t33  =  XtJ33  0n  11 


Accordingly,  X  is  a  principal  value  of  G  at  all  points  of  n,  with  the 
x3-axis  as  an  associated  principal  axis.  This  Identifies  x  as  the  transverse 
principal  stretch  on  n  and  entitles  us  to  take 

O  1 

Xj=x  on  n  .  (2.17) 

Also,  since  o33 = 0  by  the  first  of  (2.11),  one  infers  from  (2.16)  that 

o  0  ?  0  4 

p  =  Ax  -  Bx  on  n  (2.18) 


Bearing  in  mind  (2.6),  (2.9),  and  (2.16),  we  set 


1  =  tr*-GaB^  =  X1  +  A2’  J  '  detl-FaB-*  '  *1*2  * 


and  note  with  the  aid  of  (2. 3), (2. 17) ,  (2.9)  that 


(2.19) 


Note  that  x  <  1  represents  a  thinning,  and  x>  1  a  thickening,  of  the  sheet. 
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-1  0  -?  0  ?  -? 
x«  J  ^  =  1  +  J  \  I2  =  J^  +  IJ  6  . 


(2.20) 


We  now  evaluate  a  „  by  means  of  (2.7),  making  use  of  (2.18),  (2.20),  and 

ao 

thereafter  appeal  to  (2.19)  and  the  Cayley-Hamilton  theorem  in  two  dimensions 

o 

in  order  to  eliminate  G  .  from  the  resultinq  stress-deformation  relation. 

aB 

This  computation  yields 


0  000  00000^01 

-w  *BJ' K 


(2.21) 


Finally,  (2.5)  and  (2.16)  justify 


aB  ay  By 


on  n 


(2.22) 


For  convenience  we  summarize  at  this  point  in  coordinate-free  notation 
the  equations  governing  the  approximate  two-dimensional  theory  of  plane  stress 
reviewed  above.  To  this  end  we  shall  from  here  on  let  x,  y,  u  and  F,  G, 

rv  /v  /v  /v 

a,  denote  the  two-dimensional  vector  and  tensor  fields  with  the  components 

o  o  0  °  o  o 

x  ,  y  ,  u  and  F„,G„,o„,t„.  At  the  same  time  we  omit  henceforth 
a  Ja  a  aB  aB  aB  aB 

the  superscript  zero  from  symbols  that  denote  restrictions  to  n  of  functions 
originally  defined  on  the  cylindrical  region  ft. 

Thus,  (2.12)  and  (2.14)  now  become 


y(x) =  x  +  u(x)  on  n  ,  (2.23) 

iv  /v  iv  rv  iv 

diva  =  0  on  n  .  (2.24) 

Next,  keeping  (2.19),  (2.20)  in  mind,  we  introduce  the  plane-stress  elastic 
potential  U  through 

ToTT  °.t 

Fga  stands  for  the  components  of  F_  . 
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U(I,0)  =  »(!,  ,I2)  -  W(I  +  J‘2,02+  IJ'2) 

Equation  (2.21)  then  readily  leads  to 

a  =  2UTF  +  Jll.F  on  n  , 

].~  Jrw 


(2.25) 


(2.26) 


where 


F  =  Vy  ,  G  =  F  F1  , 

I  =  tr  G  =  X^  +  X^ ,  0  =  det  F  =  X^2  =  1/X  on  n  , 


(2.27) 


while  Uj  and  Uj  designate  the  corresponding  partial  derivatives  of  U(I,0). 
Further,  (2.22)  implies 


L*  2.1  on  n 


(2.28) 


Suppose  now  r  is  a  regular  arc  in  n  and  r*  =  y(r)  its  deformation 
image  in  n*  =  y(n).  Moreover,  let  n  and  n*  be  the  orienting  unit  normal 
vectors  of  r  and  r*t  respectively.  Then  (2.28)  is  easily  found  to  imply 

x  n*  *  0  on  r*  if  and  only  if  an  =  0  on  r  ,  (2.29) 

so  that  the  vanishing  of  the  nominal  tractions  on  r  is  necessary  and  suf¬ 
ficient  in  order  that  the  actual  tractions  vanish  on  r*.  Note  that  the  three- 
dimensional  counterpart  of  this  result  was  used  earlier  in  writing  the  boundary 
conditions  (2.10). 

We  also  recall  parenthetically  that  (2.26)  becomes  identical  with  the 
analogous  two-dimensional  constitutive  relation  for  plane  strain  of  a  homo¬ 


geneous  and  isotropic  compressible  material,  provided  U(I,J)  is  regarded  as 


the  plane-strain  elastic  potential. 

We  now  deduce  an  Implication  of  the  constitutive  relation  (2.26)  that 
will  play  an  essential  role  later  on.  Let  M  and  N  be  the  invariants  of 


a  o  def 1 ned  by 


Mstr(oo"^),  N  =  det  a 

'rv/v  /v 


From  (2.26),  (2.27)  follows 


a  oT  =  4U?G  +  4JU.U  ,1  + J2U^G_1 


Now  G  satisfies  the  Cayley-Hamilton  equation 


(T  -  (trG)G  +  (det  G)J_=  0  , 


so  that,  upon  post-multiplying  (2.32)  by  G  ,  one  arrives  at 


trS.‘,  =  IJ'2 


whence  (2.30),  (2.31)  yield 


M  =  (4U2  +  Uj )  I  +  81)  jUjJ 


On  the  other  hand,  post-multiplication  of  (2.26)  by  fJ  at  once  gives 


0  det  a  =  det(2UTG  +  JU,1 )  . 

rv  [r*  Jr~ 


With  the  aid  of  the  characteristic  polynomial  of  the  tensor  2UjG,  (2.35)  is 


seen  to  furnish 


Thus  M  and  N  are  analogous  to  the  scalar  invariants  I  and  J  of  the 


deformation  tensor  G 
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Jdetg,=  J2Uj  +  2JUIUJtrG  +  4Ujdetg,  ,  (2.36) 

so  that  (2.27)  and  the  second  of  (2.30)  enable  one  to  write 

N  =  2UjUjl  +  (4Uj  +  Uj)J  .  (2.37) 

Elimination  of  the  explicit  dependence  on  I  between  (2.34)  and  (2.37)  results 
in  the  identity 

2UjUj  tr  (juJ)  -  (4Uj  +  Uj)det  o  +  (4U2  -  U2)2J  =  0  ,  (2.38) 

which  is  the  desired  consequence  of  (2.26). 

We  conclude  this  section  by  considering  as  a  special  case  the  Neo-Hookean 
material,  with  which  we  shall  be  concerned  hereafter.  In  this  instance  the 
three-dimensional  elastic  potential  is  given  by 


W(I1,I2)=|(I1  -3)  , 


(2.39) 


where  y>0  is  the  material's  shear  modulus.  In  view  of  (2.25),  the  cor¬ 
responding  plane-stress  elastic  potential  takes  the  form 


U(I,J)=|(l  +  J‘2-3) 


(2.40) 


and  (2.26),  in  component  form,  thus  reduces  to 


°aBSp(FuB  FaB  =  ya,B  ’ 


X»  J 


-1 


(2.41) 


But 


JF"1  =  E  e.  F  ,  J  =  det[F  J 

Ba  ay  Bv  yv  aB 


(2.42) 
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1f  e :  is  the  two-dimensional  alternator  (ei i  =  e22  =  G12  =  "E21  =  ^‘ 

Hence  (2.41 )  yield 


aa$  U^a,6  ^  £ay£Bv^u,v^ 


(2.43) 


while  according  to  (2.28),  (2.41),  the  actual  stresses  obey 


W'Ve.v'*  s-6>  • 


(2.44) 


Substituting  from  (2.43)  Into  the  component  version  of  (2.24)  and  remembering 
the  first  of  (2.27),  one  arrives  at  the  equilibrium  equations  in  terms  of  the 
spatial  coordinates: 


7  ya  ^  ),BEoaieBv^u,v 


(2.45) 


or 


v2yl  *  (x3),ly2,2‘(x3),2y2,l  , 

v2y2S  (*3),2yl  ,r  ^3),lyl  ,2  * 


-  (2.46) 

J 


Further,  in  view  of  (2.40)  and  since  J  =  l/x,  the  relation  (2.38)  here  becomes 


X12-4(*7  +  A)N-2X6--V4M*l  =  0  ’ 

V  V 


(2.47) 


while  (2.30),  (2.43)  lead  to 


N-<let[»o6]  =  „2(A5.»-'-»\ieyaie)  , 

M"<’«8,’«8,l,2t0+>6lya,Sya.B'4x2]  ' 


(2.48) 
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The  foregoing  equations  agree  with  the  analogous  equations  In  the  paper 
by  Wong  and  Shield  [15]  J  In  particular,  (2.47)  is  in  agreement  with  Eq.(5.4) 
of  [15],  which  is  deduced  there  in  a  different  manner. 

Finally,  It  will  be  helpful  to  recall  the  response  of  a  Neo-Hookean 
material  to  a  three-dimensional  homogeneous  deformation  corresponding  to  simple 
shear  or  uni-axial  stress.  For  this  purpose  we  first  note  with  the  aid  of 
(£.3),  (2.6),  (2.8),  and(2.39)  that  the  first  of  the  constitutive  relations 
(2.7)  at  present  reduces  to 

T=uFFT-pl,  F  =  Vy  .  (2.49) 

•V  r  /*•’  IV  ^ 

A  simple  shear  parallel  to  the  plane  x3  =  0  is  represented  by  the  mapping 
y1  =x1 +  kx2,  y2  =  x2,  y3  =  x3  ,  (2.50) 


where  k  is  the  amount  of  shear  (shearing  strain).  According  to  (2.49)  this 
deformation  Induces  the  true  stresses 


12 


=  uk,  T 


23  31 


0,  t i  ^  =  yk  ,  x 


22  '  t33 


(2.51) 


provided  the  arbitrary  hydrostatic  pressure  p  is  adjusted  so  as  to  enforce 
the  vanishing  of  the  transverse  normal  stress  x33. 

The  deformation  appropriate  to  uni-axial  stress  parallel  to  the  x2-axis, 
in  turn,  takes  the  form 

y-j  *  x.|//r,  y2  *  x2,  y3=x3//T  (*>0)  ,  (2.52) 


It  should  be  noted,  however,  that  the  plane-stress  elastic  potential  employed 
in  [15]  Is  tU(I,J),  if  U  is  the  potential  defined  in  (2.40).  Further,  in 
[15]  the  equations  are  cast  in  terms  of  thickness  stress-resultants  (forces 
per  unit  length),  rather  than  in  terms  of  stresses.  Also  t  is  not  restricted 
to  be  constant  in  [15]. 


so  that  the  associated  principal  stretches  obey 


Xg  *  X,  X-|  =  X3  =  1  / 


(2.53) 


Moreover,  p  in  (2.49)  must  now  be  chosen  so  as  to  Insure  that  s  t33  =  ®’ 
In  this  manner  one  arrives  at  the  true  stresses 


T22  =  T(x)  =  u(X*"  -  X  ),  xn  =  T33  =  0  , 


c12"  t23'  t31 


L  (2.54) 


A  graph  of  the  first  of  (2.54),  depicting  the  axial  stress  t  as  a  function 
of  the  axial  stretch  X,  Is  given  In  Fig. 2. 

3.  The  global  nonlinear  Interface-crack  problem  for  two  Neo-Hookean  sheets. 
Local  analysis. 

Here  we  first  formulate,  within  the  finite  theory  of  plane  stress  for 
Neo-Hookean  materials,  the  global  interface-crack  problem  whose  local  analysis 
constitutes  our  objective.  Figure  1  is  henceforth  taken  to  refer  to  two  semi- 
infinite  Neo-Hookean  sheets,  the  interiors  of  which  occupy  the  open  half-planes 
and  H2  in  an  undeformed  configuration,  while  (x1  ,x2)  now  represent 
material  coordinates  in  this  configuration.  Further,  u-j  and  u2  at  present 
denote  the  respective  shear  moduli  of  the  two  sheets.  The  latter  are  bonded 
along  the  Xj-axls  but  for  the  crack  of  length  2s.,  which  is  assumed  to  be 
traction-free,  each  sheet  being  subjected  to  a  homogeneous  deformation  at 
infinity. 

One  is  thus  required  to  find  a  deformation  ya(x-j*x2)  that  satisfies 
the  equilibrium  equations  (2.45)  on  and  —  through  (2.43)  —  generates  a 


nominal  stress  field 
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<’aS*»k(y«.,8-l3ta/BvJ'l,,v)  °"  Hk  • 
x  =  J'1 ,  J  =  det[y  J  , 

a  f  p 

obeying  the  boundary  conditions 

oa2(Xl,0±)“°  • 

In  addition,  the  desired  deformation  is  to  conform  to  the  kinematic  loading 
conditions 

ya(xl  ,x2)  a  F^^x0 +  °{1 )  as  oo  ,  (x1,x2)£Hk  ,  (3.3) 

co(  l  \ 

in  which  P  1  are  pre-assigned  constants,  and  y  Uj.Xg)  must  meet  the 
bond  conditions 

ya(x1 ,°+) =ya(xi,0-),  aa2(x1 ,0+) = oa2(x1,°-)  ( |X| |  >  l)  .  (3.4)2 

Finally,  y  and  a  .  are  to  have  the  same  smoothness  as  u  and  o  .  In 
a  a8  a  a6 

England's  problem.3 

°o  (  1/  \ 

The  gradient  components  F^e'  of  the  homogeneous  deformation  at  in¬ 
finity  cannot  be  assigned  at  will  since  (3.3)  must  be  compatible  with  the 
bond  conditions  (3.4).  With  the  aid  of  (3.1)  one  is  thus  led  to  the  four 
conditions 

See  the  remark  following  (2.29),  according  to  which  (3.2)  assure  that  the 
deformed  crack-faces  are  free  of  actual  tractions. 

p 

The  second  of  (3.4)  is  easily  seen  to  imply  the  continuity  of  the  actual 
tractions  across  the  deformed  bonded  portions  of  the  Interface. 

See  the  regularity  assertions  following  (1.4). 


(3.1) 


(3.2)1 
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•0 )  m  p(2) 


al 


al 


''I 

>  (3.5) 


j 


where 


s  = 


det[. 


>(k) 

aB 


]>0 


( 3 . 6 ) 1 


The  relations  (3.5)  are  a  counterpart  in  the  present  nonlinear  settinc 
the  loading  constraint  (1.15). 

If  3  u2  (s*l),  so  that  one  is  dealing  with  a  crack  in  a  single 
all-aroung  homogeneous  Neo-Hookean  sheet,  (3.5)  are  automatically  satisfied, 
provided 


?0)  =  ?(2)  =  c 

aB  aB  '  aB 


(3.7) 


Moreover,  the  particular  assignment 


oo 

X2  >  1 


(3.8) 


00  00 


corresponds  to  a  Mode  I  loading  with  the  principal  stretches  ,  x^  at 
infinity,  while 


[y- 


1  k 
0  1 


(3.9) 


constitutes  a  Mode  II  loading,  k  being  the  amount  of  shear  applied  at 


T 


The  inequality  in  (3.6)  follows  from  the  inequality  in  (2.41). 
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[  infinity. 

With  a  view  toward  an  asymptotic  analysis  of  the  elastostatic  field 
near  the  right-hand  crack-tip,  it  is  expedient  to  regard  the  spatial  coordi¬ 
nates  y^  from  here  on  as  functions  of  the  local  (material)  polar  coordinates 
(r,e)  (see  Fig . 1 )  and  to  introduce  the  complex  spatial  coordinate 

y(r,e)  =  y1(r,e)  +  iy2(r, s)  (r  >  0,  -it  se  in)  .  (3.10) 


One  then  has 


J  =  detC4,8]ar"1lin{lf  le}  *  X  =  J  1  *  (3J1) 

provided  3y/3r  stands  for  the  complex  conjugate  of  3y/3r.  The  equilibrium 
equations  (2.45)  are  equivalent  to 


b!. 


V  y 


|^r  39 


36  3r 


3r^  r  3r  rc  30^ 


y  (3.12) 


while  the  constitutive  relations  (2.43)  lead  to 


all  +  ia21  "  u[fr^cos  9  +  ix3s,in  ~  f  fe  (sin  6  -  ix3cos  e)j  , 

a12  +  io22  =  u[lr  8  ‘  ix3cos  9)  +  ||- (cos  6  +  ix3sin  e)] 


(3.13) 


From  (3.13)  one  deduces  formulas  for  the  polar  components  of  nominal  stress: 
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0rr+1o8r 


*  ue~i0-|  +  U3  , 

[ Sr  r  36  J  ’ 


cr6  +  1 


-ie  f  1  ay  .,3  ay\ 
laee  =  ue  \r  96  '  u  a?)  • 


A 


(3.14 


,1 


Further,  (2.48)  at  once  yield 


>  (3.15) 


N  ■  <tettoa#] -u2{xi-x-1  -  x3  (l|^|2-  ^-llil2)}  • 


H  =  %6°»e"“2{(’**6>(1^1  +7I“I  )"41' 


Equations  (3.12)  hold  for  0<8<ir,  as  well  as  for  -Tt<e^0.  On  the 
other  hand,  (3.13),  (3.14),  and  (3.15)  apply  to  0&6 &tt  with  u  -  and 
to  -it  40 sO  with  u=p2*  Moreover,  because  of  (3.2),  (3.13),  one  has  the 
boundary  conditions 


fl  -U3  &] 
Lr  ae  arj 


(r,±Ti) 


=  0  (0  <  r  <  2l)  , 


(3.16) 


whereas  (3.4),  (3.13)  furnish  the  bond  conditions 


y(r,0+)  =y(r,0-)  , 


s  [1 9X  .  ix3  ax]  =[liZ_n3  3*.l 
Lr  ae  3rJ(r>0+)  lr  30  3r  J (r-.Q-J 


(r  >  0)  , 


(3.17) 


and  analogous  bond  conditions  at  e  =  ±ir. 

Taking  for  granted  the  existence  of  a  solution  to  the  global  interface- 


The  subscripts  in  the  left  members  of  (3.14)  are  exempted  from  the  range 
and  summation  conventions  adopted  earlier  for  Greek  indices. 
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crack  problem  stated  at  the  beginning  of  this  section,  we  now  assume  that 
this  solution  admits  an  asymptotic  representation  of  the  form 

y{r  ,0)  ~i  +  rmv(e)  as  r-+0  (Osqstt)  ,  'l 


v(r.al  +  r 


.m 


(3.18) 


1 


m  *  a  +  ie,  m'*  a'  +  ie'  (a, a'  and  8,8'  real)  ,  (3.19) 

so  that  the  exponents  m  and  m'  are  permitted  to  be  complex  numbers, 
while  v  is  a  complex-valued  function,  twice  continuously  differentiable  on 
[0,ir]  and  on  [-*,0],  that  fails  to  vanish  Identically  on  either  of  these 
two  intervals.  Further,  we  stipulate  that 


a  =  Rem  >  0,  ex'-  Re  m'>  0 


(3.20) 


accordingly  the  displacements  are  required  to  be  continuous  as  r+0  and  the 
elastostatlc  field  at  hand  Is  taken  to  have  been  normalized  by  means  of  a 
suitable  rigid  translation  so  as  to  keep  the  crack-tip  at  r  =  0  fixed. 

Note  that  (3.18),  (3.19),  (3.10)  furnish 


y^r.e)  ~i  +  ra(v1  (e)cos(e  log  r)  -  Vgtejsin^  log  r)J 
~<"a{v^(e)sin(B  log  r)  +  Vg^costs  log  r)] 


i  (3.21 )2 


when  OiQiir;  if  -tsq&O,  a  and  b  in  (3.21)  are  to  be  replaced  by  a 


Here  and  In  the  sequel  the  asymptotic  equality  symbol  Is  used  in  the 
following  connotation:  the  first  of  (3.18)  means  that  for  O^esit, 

y(r,e)  =  t  + rmv(e) +  o(rm)  as  r-*-0  , 

where  o(rm) * o( |rm| ) » o(ra) ,  If  a  =  Rem. 

2 

All  asymptotic  ..-qualities  and  order-of-magnitude  estimates  are  henceforth 
understood  to  refer  to  the  limit  as  r-<-0. 
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and  b',  respectively.  Equations  (3.21)  exhibit  the  structure  of  the  oscil¬ 
latory  behavior  admitted  by  the  Ansatz  (3.18).  Clearly,  the  right-hand 
members  in  (3.21)  are  non-oscillatory  if  and  only  if  m  and  m' .  are  both 
real  (b=b'3Q).  Observe  also  that  (3.18)  encompass  the  possibility  that 
y.|(r,8)-Jl  and  y2(r,e)  are  of  different  dominant  orders  of  mangitude  as 
r+0,  for  a  fixed  value  of  9.  For  example,  according  to  (3.21),  a  =  0 
together  with  v2  3  0  on  [0,tt]  gives  y-j  (r,e)  -  I ~rav-|  (e),  but  y2(r,e) 
in  this  instance  vanishes  to  a  higher  order  as  r->-Q. 

To  (3.18)  we  adjoin  an  additional  asymptotic  assumption  concerning  the 
transverse  stretch  x  by  requiring  that 

X(r,e)  =  1/J(r,9)  =  0(rq)  as  r+0  (-ir&e&ir)  for  some  q>0  .  (3.22) 

At  the  same  time  we  demand  the  validity  of  the  asymptotic  equalities  resulting 
from  two  consecutive  formal  differentiations  of  (3.18)  with  respect  to  r  or 
9,  as  well  as  the  "differentiability"  of  (3.22). 1  Equation  (3.22)  implies 
that  J(r,e)  becomes  unbounded  as  r-*-0  and  hence  anticipates  that  not  all 
deformation-gradient  components  remain  bounded  in  this  limit,  so  that  «<1, 
ot'<l.  Since  x(r,e)  tends  to  zero  as  r-<-0,  in  view  of  (3.22),  the  deforma¬ 
tion  prevailing  in  the  vicinity  of  the  crack-tip  entails  an  extreme  thinning 
of  each  sheet.  Evidently  (3.22)  excludes  all  regular  deformations  -  in  parti¬ 
cular  homogeneous  deformations  -  from  our  present  considerations. 

Aiming  at  the  lowest-order  asymptotic  near-field  behavior,  we  now  seek 
to  determine  v(r,e),  as  well  as  the  exponents  m,  m'  with  the  smallest  real 
parts,  consistent  with  (3.18),  (3.20),  (3.22),  the  governing  equilibrium 
equation  (3.12)  and  the  accompanying  boundary  and  bond  conditions  (3.16), 

Thus  vx(r,9)  =  0(rq_1 )  s  r-*-0  (0<je|<»). 
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(3.17).  On  account  of  (3.18)  and  (3.22),  the  nonlinear  x^-terms  in  (3.12), 
(3.16),  (3.17)  do  not  come  into  play  in  the  lowest-order  analysis.  Upon 
entering  (3.12),  (3.16)  with  (3.18),  (3.22),  one  readily  arrives  at 

v  +  m^v*0  on  (0,ir),  v  +  (m'  )^v  *  0  on  (-n,0)  ,  (3.23)1 

v(it)  a  v(-ir)  *  0  .  (3.24) 


From  (3.23)  follows 


v(e)  =  a  sin  me  +  b  cos  me  (0<e*ir)  , 

v(e)  =  a'sinm'e +  b'cosm'e  (-ir£8<0)  , 


V  (3.25) 


In  which  a,  b  and  a',  b'  are  complex  constants.  On  the  other  hand,  (3.18), 
(3.22)  in  conjunction  with  the  bond  conditions  (3.17)  and  (3.25)  give 


s  ma 


br!TI  +  o(ra)  »b'  r"1  +o(ra)  , 
rm_1  +o(ra-1)=mVrm''1  +  o(ra'-'1) 


l  (3.26) 


If  at*',  (3.26),  (3.20)  imply  that  a,  a',  b,  b'  vanish  simultaneously, 
which  is  Inadmissible  since  in  this  instance  v  =  0  on  [-*,*].  Therefore, 


a  =  Re  m  *  a*  =  Re  m/ 


(3.27) 


Further,  (3.26),  (3.27)  necessitate  that 


^Here  and  In  what  follows  superior  dots  connote  differentiation  with  respect 
to  the  polar  angle  e. 
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b  +  o(l)  =  b'ri(e-6/)  , 

s  m  a  +  o(l )  =  m'a'r1'  ^  , 


(3.28) 


and  the  left  members  in  (3.28)  tend  to  b  and  sma,  respectively,  as  r->-0, 
whereas  the  right  members  fail  to  possess  limits  unless  b=b'.  Consequently, 
one  draws  from  (3.27)  and  (3.28)  that 


m  =  m' ,  s  a  =  a' ,  b  =  b' 


(3.29) 


Combining  (3.29)  with  (3.25),  invoking  (3.24),  and  recalling  that  m=0  is 
precluded  by  (3.20),  one  has 


a  cos  mn  -  b  sinirnr  =  0  , 

s  a  cos  mrr  +  b  sinm*  *  0 


(3.30) 


Thus,  since  s>0  and  a  =  b=0  is  ruled  out  by  v^O  on  [0,ir], 


sin2mir  =  0  , 


(3.31) 


whence  m  is  real .  The  smallest  positive  root  of  (3.31)  is 

m  =  a  =  ^  >  (3.32) 


and  (3.30),  (3.32)  lead  to  ba0,  while  the  constant  a  remains  arbitrary. 
Also,  (3.32),  (3.29),  (3.25),  and  (3.18)  yield 


y(r,e)  +  r1/,2v(e)  (-nsesn)  , 

v(e)  -  ah(e)sin|-  (-ns.es.Tt)  , 


■'i 


(3.33) 


in  which  a t  0  Is  an  arbitrary  complex  constant  and  h  is  the  step-function 
defined  by 


h ( 0 )  =  1  (OsesiTr),  h(e)  =  s  {-irie<0)J  s  =  y1/u2  .  (3.34) 


On  setting  a  =  a-|  +  ia2  and  referring  to  (3.10),  one  finally  concludes  that 

y1  (r,e) +  a1r1/2h(e)sin  |  , 

y2(r,e)  ~a2rly/2h(e)sin  |  (-irsiasir)  , 


>  (3.35) 


2  2 

where  ,  a2  are  arbitrary  real  constants  with  a1+a2>0. 

One  confirms  with  the  aid  of  (3.11)  that  the  Jacobian  determinant  of  the 
right  members  in  (3.35)  vanishes  identically.  Therefore  the  lowest-order 
approximation  to  the  local  deformation  (near  the  crack-tip)  obtained  above 
does  not  constitute  a  mapping  that  Is  one-to-one.  Furthermore,  since 
X  * 1/J ,  the  one-term  asymptotic  solution  (3.35)  fails  to  supply  an  estimate 
for  the  transverse  stretch.  These  inadequacies  make  it  necessary  to  seek  at 
least  a  two- term  approximation  to  the  elastostatic  field  in  the  vicinity  of 
the  crack-tip. 

In  preparation  for  this  task  we  show  first  that  (3.33),  together  with 
(2.47)  and  (3.22),  in  fact  enable  one  to  deduce  an  estimate  for  x(r,±w)  as 
r-*-0,  i.e.  for  the  transverse  stretch  along  the  faces  of  the  crack.1  Indeed, 
the  boundary  conditions  (3.2),  referred  to  the  polar  coordinates  (r,e),  give 

N(r,±ir)  =  det[a  (r,±ir)]  =  0  (0<r<2t)  ,  (3.36) 

dp 

This  appeal  to  the  identity  (2.47)  was  suggested  by  Wong  and  Shield  [15],  who 
use  (2.47)  in  a  similar  manner  in  their  treatment  of  a  Mode  I  problem  for  a 
crack  in  an  all-around  homogeneous  Neo-Hookean  sheet. 
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so  that  (2.47)  supplies 

^  =  1+X12-2X6  at  8  =  -n  (0  <  r  <  2s.)  .  (3.37) 

P1 


Next,  (3.33),  (3.22),  and  the  second  of  (3.15)  are  found  to  imply 


M  =  0  .0 
ad  06 


4r 


(-nsesii),  a  f  0 


(3.38) 


Keeping  in  mind  that  X(r,@)  =  o(1)  for  -tsssi  and  combining  (3.38)  with 
(3.37),  as  well  as  with  its  companion  for  0  =  -it,  one  readily  confirms  that 


X(r,ir) 


x(r,-ir) 


as 


r  -*■  0  , 


(3.39) 


which  is  consistent  with  (3.22)  and  necessitates  qsl/4.  Guided  by  (3.39), 
we  now  refine  the  original  a  priori  assumption  (3.22)  by  anticipating  that 


x(r ,e)  *  0(r^  )  as  r -*•  0 
vx(r,e)  =  0(r'1,/4)  as  r-*-0 
It  is  helpful  to  observe  that  (3.40),  (3.33) 

X3|^  =  °(r1/4)  as  r-*-  0 


{-ir  s  0  Sir)  , 


(0<  |  0  |  <  it) 
imply 

(-ir  s  0  Sir)  , 


>  (3.40) 


j 


(3.41) 


i 

r 


/ ay  ax3  ay  ax3N\ 

\3r  80  00  8r  / 


=  0(r-3/4) 


as  r+  0 


(-ir  S  0  Sir) 


(3.42) 


Proceeding  to  the  second-order  asymptotic  analysis,  we  set 
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y(r ,e)  +  r^2v(e)  +  rn  w(e)  (O&esu)  , 

y(r,e)  ~t  +  r1/2v(e)  +  rn  w(e)  (-nsesO) 


(3.43) 


Here  v  is  known  from  (3.33),  n  and  n'  are  as  yet  unknown  complex  ex¬ 
ponents  obeying 


Y  =  Re  n  >  -^  ,  y'  =  Re  n  '>  ^  » 


(3.44) 


and  w  is  an  initially  undetermined  complex-valued  function  that  is  not 
permitted  to  vanish  identically  on  [0 ,-rr]  or  [-ir,0].  Further,  w  and  (3.43) 
are  subject  to  the  same  differentiability  hypotheses  made  previously  in  con¬ 


nection  with  v  and  (3.18). 

Our  current  objective  is  to  find  w  on  [-■*,*]  and  n,  n' ,  with  the 
smallest  real  parts  admitted  by  (3.44),  consistent  with  (3.12),  (3.16),  (3.17), 
and  (3.40).  In  view  of  (3.41),  (3.42),  (3.43),  and  because  v  satisfies 
(3.23),  (3.24),  the  equilibrium  equation  (3.12)  and  the  boundary  conditions 
(3.16)  require 


rn-2[«(e)  +  n2w(e)]  +  o(rn'2)  =  0(r'3/4)  (0<e<*)  , 

rn' 2[w(0)  +  (n')2w(e)]  +  o(rn' 2)  =0(r'3/4)  (-*  <  0  <  0) 


-  (3.45) 


If 


rn‘1;(^)  +  o(rn-1)=0(r1/4)  , 

rn ’^w(-ir)  +  o(rn  "1)=0(r1^4) 

Re  n  <  | ,  Re  n'  <  |  , 


|  (3.46) 


(3.47) 
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as  we  shall  temporarily  take  for  granted,  (3.45),  (3.46)  evidently  imply 
(3.23),  (3.24)  -  and  hence  (3.25)  -  with  v,  m,  m'  replaced  by  w,  n,  n' , 
respectively.  Also,  (3.47)  injure  that  the  nonlinear  x  -terms  in  the  bond 
conditions  (3.17)  likewise  fail  to  enter  the  second-order  analysis  explicitly. 

An  argument  strictly  parallel  to  that  following  (3.25)  easily  completes 
the  second-order  asymptotic  solution.  In  particular  one  finds  in  this  manner 
that 

n  =  n' ,  sin  2n*  =  0  ,  (3.48) 

whence  n  and  n'  are  real .  The  smallest  root  of  (3.48)  conforming  to  (3.44) 
is 

n  =  n'=  y  =  y1  =  1  ,  (3.49) 


and  since  (3.49)  is  consistent  with  (3.47),  there  is  no  need  to  explore  the 
complementary  ranges  yis5/4,  y's5/4.  The  second-order  approximation  thus 
emerging  is  given  by 


y(r,9)  +  r1/,2v(e)  +  rw(e)  (-tts9s.it)  , 

v(e)  =  a  h(e)sin  ,  w(e)  =  b  cos  e  (-use  sir)  , 


(3.50) 

J 


In  which  a  f  0  and  bf*0  are  otherwise  arbitrary  complex  constants,  while 
h  again  stands  for  the  step-function  Introduced  in  (3.34).  From  (3.50), 
(3.10)  finally  follows: 


1/2  0 

y-j  (r,8)  +  a^r  '  h(e)sin ^  +  ^r  cos  e 


y2(r,9)  ~a2r1/2h(e)sin!  +  b2rcos  6  (-irsesn)  , 


(3.51) 
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with 


7  7  7  7 

a,  +  a,  >  0,  bf  +  b~>0 


(3.52) 


ap  and  bi«  b2  be^n9  otherwise  arbitrary  real  constants. 

Upon  using  (3.11)  to  compute  the  Jacobian  determinant  on  the  basis  of 
the  two -term  asymptotic  solution  (3.50),  one  arrives  at  the  estimate 

J(r,e)  ~Ar  ^^(ejcosj  (-irs.9&ir)  ,  (3.53) 

in  which 

A*  ^(a^k]  “  a-^)  *0  ,  (3. 54 ) 

the  inequality  being  a  consequence  of  the  requirement  J>0.  The  approximation 
(3.53)  gives  merely  J(r,e)  *  ofr"1^2)  (-usosn)  unless  A>0,  as  we  shall 
assume  until  further  notice.1  From  (3.53)  one  draws  the  asymptotic  charac¬ 
terization  of  the  transverse  stretch  x  =  l/J: 

I/? 

x(r’9)  ~AH{eJc'os(e/2)  (_1t<0<1t)  lf  A>0  *  (3*55) 

which  is  consistent  with  the  a  priori  assumption  (3.40)  employed  in  the  deri¬ 
vation  of  (3.50).  Equations  (3.53),  (3.55),  (3,34)  reflect  the  discontinuity 
in  J(r,e)  and  x(r,e)  at  9  =  0  if  s  f  1  (u-j  ^  vig ) * 

Although  (3.55)  and  (3.39)  together  furnish  estimates  of  x(r,9),  as 
r  +  Q,  for  each  fixed  9  in  the  complete  range  [-ir.n],  the  approximation 
given  by  (3.55)  evidently  deteriorates  severely  near  the  crack-faces  6  =  ±ir. 

In  order  to  obtain  a  crack-tip  estimate  for  *  that  is  free  of  this  deficiency 

The  special  case  A=0  will  be  treated  in  the  next  section. 
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and  that  at  the  same  time  reveals  the  nature  of  the  transition  from  (3.55) 
to  (3.39)  as  0-*±ir,  one  may  once  again  appeal  to  the  Identity  (2.47)  used 
earlier  to  deduce  (3.39).  Since  (3.22)  Implies  x(r.e)-o(1)  as  r-0, 

(2.47)  near  r«Q  Is  approximated  by 

M(r,e)x4  +  N(r,@)x  -  u2  ■  o  (Oiesn)  .  (3.56)1 


This  suggests  that 

A(r,e)  ~xo(r,e)  as  r-*-0  (Osesir)  ,  (3.57) 

where  \Q  is  the  unique  positive  root  for  x  of  (3.56),  as  is  not  difficult 
to  confirm  rigorously.  On  setting 

x-i^3/2M’1/4N,  n»u^1/2M1/4x  ,  (3.58) 

(3.56)  becomes 

n4  +  xn  -  1  =  0  .  (3-59) 

A 

Let  na<p(x)  be  the  unique  positive  root  of  (3.59),  so  that  (3.58),  (3.57) 
give 

x(r,e)~u]/2M'1/49(y^3/2M"1/4N)  (0*0**)  .  (3.60) 

The  desired  estimate  for  x  Is  obtained  by  inserting  small -r  approxi¬ 
mations  for  M(r,e)  and  N(r,e)  Into  (3.60).  Such  an  approximation  to 

^eTketch  in  detail  only  the  analysis  pertaining  to  the  upper  sheet  (0*e*ir). 
Recall  from  (2.48)  that  M-  aa0oaB >  °*  N»det[oa6]. 

ZAn  explicit  elementary  representation  for  cp(x)  is  readily  established,  but 
serves  no  particular  purpose  here. 
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M(r,e)  Is  available  from  (3.38).  On  the  other  hand,  an  analogous  estimate 
for  N(r,e)«det[o  J  Is  computable,  without  recourse  to  (3.55),  from  (3.50) 
and  (3.13)  by  relying  merely  on  (3.40).  In  this  manner  one  obtains 

N(r,e)~u2Ar'1/2cos|  (0  * e  <  *),  N(r,n)  *  o(r’1/2)  .  (3.61 


Substitution  from  (3.38),  (3.61)  into  (3.60),  along  with  a  parallel  argument 
applied  to  the  lower  sheet,  results  in 


x(r,e)  ~r1/,4g(e)tp(r*1/,4Ag(e)cos(e/2))  , 

g(e)  *  J2/  |a)  h(e)  (-ime**),  A>0  . 


(3.62) 


j 


Under  suitable  assumptions  concerning  the  uniformity  in  e  of  the  Ansatz 
(3.43),  one  can  show  that  the  estimate  (3.62),  the  validity  of  which  is  con¬ 
tingent  upon  our  current  supposition  that  A>0,  holds  uniformly  for 

-IT  £  8  S>  IT  , 

Finally,  we  note  that  one  recovers  from  (3.62)  through  appropriate 
specialization,  the  previous  estimates  (3.39)  and  (3.55),  which  apply  to 
e*±ir  and  separately.  To  see  this  we  infer  from  (3.59)  that 


<p(0)  *  1 ,  <p'  (x)  <  0  (-00  <  x  <  00)  ,  «p(x)  ~  -  as  x-*-®  -  (3.63) 

The  first  of  (3.63)  reveals  that  (3.62)  reduces  to  (3.39)  for  93±ir;  the 

last  of  (3.63)  leads  immediately  from  (3.62)  to  (3.55). 

Although  we  have  not  determined  in  detail  any  further  terms  in  the  ex- 

5/4 

pension  (3.50),  we  have  ascertained  that  additional  terms  of  orders  r  , 
Here  A  is  again  the  auxiliary  amplitude  parameter  defined  in  (3.54). 
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r37^,  and  r774  can  be  obtained  by  arguments  parallel  to  those  used  in  this 
and  the  preceding  section.  If  one  attempts  to  continue  the  expansion  beyond 
the  term  of  order  r77^,  however,  one  encounters  difficulties  since  at  that 
stage  the  x '‘-terms  in  the  equilibrium  equation  (3.12)  begin  to  enter  ex¬ 
plicitly  into  the  asymptotic  analysis.  Indeed,  it  appears  from  the  structure 
of  the  uniform  approximation  (3.62)  for  X  that  such  higher-order  terms  can 
no  longer  be  of  the  "factored"  form  r  f(e). 

4.  Rotation  of  the  deformed  composite  sheet.  The  degenerate  case  A=0. 

The  complex  version  (3.50)  of  the  two-term  asymptotic  approximation  de¬ 
duced  in  Section  3  involves  the  complex  amplitude  parameters  a  and  b 
(a?4  0,  b><  0),  the  determination  of  which  eludes  the  local  analysis  carried 
out  there.  The  values  of  these  two  parameters  are  bound  to  depend  on  the 
material  parameters  u1  and  the  crack-length  2e,  as  well  as  on  the 
particular  loading  at  infinity,  which  has  not  entered  the  derivation  of 
(3.50).  Indeed,  these  local  results  are  not  confined  to  semi-infinite  sheets 
and  encompass  essentially  arbitrary  loadings  applied  to  the  outer  boundary. 

We  show  presently  that  (3.50)  may,  without  loss  of  generality,  be  simplified 
by  subjecting  the  deformed  composite  sheet  to  a  rigid  rotation  about  the 
crack-tip  at  r  =  0.  Moreover,  this  transformation  -  as  will  become  evident  — 
renders  the  intrinsic  content  of  the  asymptotic  solution  (3.50)  more  trans¬ 
parent. 

If  y(r,e)  is  the  complex  spatial  coordinate  (in  the  original  fixed 
Cartesian  frame)  of  a  material  point  (r,e)  after  such  a  rotation  through 
an  angle  «J,  one  has 


y(r,e)  =  n  +  e^y (r,e)  -  i} 


(4.1) 
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we  now  take 


where  e  1$  the  argument  of  the  complex  parameter  a,  whence 

a  =  a1 +  ia2=  |a|e1e  . 

For  this  choice  of  9  ,  equations  (4.1),  (3.50)  imply 


y(r,e)  +  §r^h(e)sin  +  br  cos  6  (-irsesn)  , 
a  =  i  |  a  | ,  b=iab/|a|  (a^o,  b^O)  . 


From  (4.5)  follows 


**«1  +  1»2*  *1  =0’  S2=  >0  * 
b*b] +1b2,  b]  *2A/|a|iO,  b2  =  2B/|a|  , 

where 

A  =  ^(a2bi -aib2)  i°,  B-^(a1az  +  b^)  , 


"1 


(4.2) 


(J.3) 


(4.4) 

(4.5) 


(4.5) 


(4.7) 


so  that  A  is  the  non-negative  auxiliary  amplitude  parameter  originally 
Introduced  in  (3.54).  Also,  (4.6)  give 


1*1  =  1*1 .  A  =  ^(a2b1  -  a^g)  «  A 
I b |  ■  |b| ,  B  =  ^(a1b1  +a2B2)  =  B 


■\ 

► 


j 


(4.8) 


which  assert  the  invariance  of  |a|,  |b|.  A,  and  B  under  the  rotation  (4.1). 
To  avoid  cumbersome  notation  we  henceforth  write  y,  a,  b  in  place  of 
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y,  a,  b  and  thus  note  that  (4.4),  on  account  of  (4.6),  leads  to 


y1  (r,e)  ~i  +  b^rcos  e  , 

1/2  6 

y2 ( r » 9 )  ~a2r  h(a)sin  j  +  b2r  cos  a  (-iseiir) 


(4.9) 


where 

a2>0,  b^O  (b.,  =  2A/a2)  ,  (4.10) 

and  &2>  b-j ,  t>2  are  otherwise  arbitrary  real  constants.  As  is  now  apparent, 
i.o  generality  is  lost  in  the  original  two-term  asymptotic  approximation  (3.51) 
by  setting  a-|=  0.  while  subjecting  a2  and  b-j  to  the  inequalities  in  (4,10) 
In  view  of  (4.10),  the  exceptional  case  A=0,  which  was  excluded  in 
establishing  the  estimates  (3.55)  and  (3.62)  for  the  transverse  stretch  x, 
corresponds  to  b^  =  0  in  (4.9).  On  the  other  hand,  when  b-j  =  0,  the  approxi¬ 
mating  deformation  supplied  by  (4.9)  is  no  longer  one-to-one1  and  the  current 
approximation  to  y-|(r,e)  degenerates  into  the  weak  estimate 

y-j  (r  ,e)  -  z  +  o(r)  as  r->-0  (-ttssstt)  ,  (4.11) 

which  is  inadequate.  The  case  A=0  (b^  =0)  thus  entails  an  essential  de¬ 
terioration  of  the  asymptotic  solution  (4.9)  and  necessitates  higher-order 
considerations. 

Motivated  by  the  foregoing  observations  and  bearing  in  mind  (4.9),  (4.10) 
with  b-j  =0,  we  now  adopt  the  three-term  asymptotic  Ansatz: 


^Slnce  0  Is  Invariant  under  the  rotation  (4.1),  the  Jacobian  determinant  of 
this  mapping  Is  given  by  the  right  member  in  (3.53),  which  vanishes  for  A=0. 


-41- 


y(r,e)  +  la^^sin^- +  ib2r  cos  0  +  rpz(e)  (OsesTr)  , 
y(r,e)  ~i  +  isa2r1/,2sin|-  +  ib2r  cos  e  +  rp  z(e)  (-tsesO) 


-\ 


l  (4.12) 


’J 


where 


a2  >  0.  Re  p  >  1 ,  Re  p'  >  1  , 


(4.13) 


while  the  as  yet  unkiown  complex- valued  function  z  =  z^  +  iz2  must  not  vanish 
identically  on  [0,ir]  or  C-tt,0]  and  is  required  to  be  twice  continuously 
differentiable  on  each  of  these  two  intervals.^ 

From  (4.12)  one  gathers 


v2y  ~rP“2(z  +  p2z}  (0  <  ©  <  r)  , 
v2y~rp  _2{z  +  (p' )2z}  (-tt<6<0)  . 


l  (4-14) 


Consequently,  the  equilibrium  equation  (3.12)  together  with  the  estimate 
(3.42)  enables  one  to  conclude  that 

z  +  p2z  =  0  (O<0<tt)  if  1  <  Re  p  <  5/4  ,  (4.15) 

z+(p')2z  =  0  (-T-  <  6  <  0)  if  1  <  Re  p'  <  5/4  .  (4.16) 

Further,  on  applying  the  boundary  conditions  (3.16)  to  (4.12)  and  taking 
account  of  the  estimate  (3.41),  one  has 

z(ir)«0  If  1  <  Re  p  <  5/4  ,  (4.17) 

^e  take  for  granted  the  "differentiability"  of  the  asymptotic  identities 
(4.12). 
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£(-*)  =  0  if  1  <  Re  p'  <  5/4  .  (4.18) 

Next,  imposing  the  transition  conditions  (3.17)  on  (4.12)  and  appealing  once 
again  to  (3.41),  we  see  that 

z(0+)  =  z(0+)  =  0  if  1<  Rep  <5/4,  Rep<Rep'  ,  (4.19) 

z(O-)  =  z(O-)  =  0  if  1  <  Rep'<  5/4,  Rep'<  Rep  .  (4.20) 

But,  according  to  (4.15),  (4.17),  (4.19),  the  inequalities  in  (4.19)  demand 
that  z  =  0  on  [0,ir],  which  is  inadmissible;  similarly,  the  inequalities  in 
(4.20),  in  view  of  (4.16),  (4.18),  (4.20),  necessitate  z=0  on  [-tt,0], 
which  is  also  unacceptable.  On  the  other  hand,  the  transition  conditions 
(3.17)  reduce  to 

z(0+)  =  z(0-),  sz(0+)  =  z(0-)  if  1  <  Re  p  -  Re  p  <  5/4  ,  (4.21) 

and  in  this  event  (4.15)  to  (4.21)  require  that  z  =  0  on  [-tt.tt]. 

The  preceding  conclusions  entitle  us  to  claim  that  Re  pi 5/4,  Rep'i5/4. 
Moreover,  since  our  present  objective  is  an  asymptotically  consistent  three- 
term  approximation  of  the  form  (4.12)  in  which  the  exponents  p  and  p' 
have  the  smallest  real  parts  admitted  by  (4.13),  there  is  no  need  to  consider 
Rep >5/4,  Rep'>5/4  If  we  succeed  in  determining  z(e)  on  the  supposition 
that  Re  p  =  Re  p' =  5/4.  We  are  thus  led  to  set 

p  =  |+ 15,  p'  =  |+ 16'  (6,6'  real)  .  (4.22) 

Calculating  the  Jacobian  determinant  J(r,e)  from  (4.12)  by  means  of 
(3.11)  and  using  (4.22),  one  arrives  at 


_ _ 
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provided 


u(r,e) 


(0  <  9  S it)  , 


(4.23) 


o)(r,e)  =  Re{(|- -  i5)r"i6z(0)cos  |  -  rl6z(9)sin|-}  .  (4.24) 

If  S?0,  Eq$.(4.23),  (4.24)  are  easily  found  to  yield  two  alternative  impli¬ 
cations:  either  J(r,e)  changes  sign  infinitely  often  as  r-+0  for  each  e 
in  (0,it]  o£  z(9)  is  such  that  w(r,9)  vanishes  identically  on  this  interval. 
The  first  of  these  eventualities  is  precluded  by  the  requirement  J(r,e)>0; 
the  second  alternative  leads  to  r^J(r,e) 0  as  r-*-0  and  hence  gives 
r“^x(r,e) oo  as  r-*-0,  for  each  8  in  (0,it],  which  is  incompatible  with 
(3.40).  Therefore  6  =  0  and  p  must  be  real.  One  shows  analogously  that 
p'  has  to  be  real  as  well,  so  that  (4.22)  now  reduce  to 

p  =  p'=  5/4,  6=  6'=0  .  (4.25) 

With  (4.25)  in  force  one  infers  from  (4.23),  (4.24)  and  their  counterpart  for 
-tt  s  a  <  0  that 

J(r,e)  ~r’^ui(e)  (-irses*)  ,  (4.26) 

w(e)  =  ^a2h(8)(|z1(e)cos|  -  z1(e)sin|)  ,  (4.27)1 

in  which  h(e)  is  the  step-function  defined  In  (3.34)  and  z^(e)  is  the 
real  part  of  z(e).  Since  x(r,e) =  1/J(r,e) ,  a  zero  of  w(e)  anywhere  on 
[ — tt , it ]  would  contradict  (3.40),  so  that 

As  a)(r,e)  is  now  independent  of  r,  we  write  u>(e)  in  place  of  w(r,0). 
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A(r,e)  ~r  '  V<4e) ,  u{e)>0  (-irsesir)  . 


(4.28) 


We  are  now  in  a  position  to  use  the  equilibrium  equation  (3.12)  in  con¬ 
junction  with  (4.14),  (4.25).  (4.27)  to  obtain  the  differential  equation 
governing  2(6).  In  contrast  to  the  lower-order  analyses  carried  out  in 
Section  3,  the  right  and  left  sides  of  (3.12)  now  balance  asymptotically, 
and  one  arrives  at 


z  +  yg-z= — ^  h ( 0 ) { jjp cos  ^  ~  “  sin  |-}  ( 0  <  j  0 1  <  -n- )  , 

2(o 


(4.29) 


where  w  is  related  to  z  through  (4.27).  The  boundary  conditions  (3.16), 
in  turn,  lead  from  (4.12),  by  virtue  of  (4.25)  and  (4.28),  to  the  requirements 


I  *5  1  3 

z(-ir)  a  2"Sa2°J  (“*)»  z(n )  =  -  £ a^u)  (n) 


(4.30) 


while  the  transition  conditions  (3.17)  are  found  to  yield 


z(0+) =  z(O-) ,  sz(0+) =  z(O-)  . 


(4.31) 


Since  the  right  members  in  (4.29)  and  (4.30)  are  real-valued  functions. 
It  follows  with  the  aid  of  (4.31)  that  z2 =  Imz  must  satisfy 


h*uz2s°  *0<  <ir)  » 


(4.32) 


subject  to  the  boundary  and  transition  conditions 


Z2(-ir)  *  Z2(ir)  =  0,  Z2(0+)  =  Zz(0-),  SZ2(0+)  =  ZgfO- ) 


(4.33) 


This  homogeneous  liner  boundary-value  problem  admits  only  the  trivial  solution 


“vis .... 


The  boundary-value  problem  governing  z1  s  Re  z  is  obtained  by  taking 
real  parts  in  (4.29),  (4.30),  (4.31).  The  resulting  problem  can  be  cast 
into  a  more  convenient  form  by  means  of  a  suitable  rescaling  of  z-j .  To 
this  end  we  set 


2]  (e)  =  23/231/4a‘1/2t(e)  (-*  se  s„)  (4.35) 

and  note  first  that  (4.27)  now  becomes 

u>(e)  =  2"3/231/4aJ/2h(9)«j.(e)  (-^e^)  ,  (4.36) 

provided 

^(e)  =  5c(e)cos I  -  4c(e)sin|-  (-irsesn)  .  (4.37) 

Further,  in  view  of  (4.36)  and  the  positivity  of  u>(e),  one  has  to  stipulate 
that 

ii>(e)>0  (-TrsesTi)  .  (4.38) 

On  equating  real  parts  in  (4.29)  and  making  use  of  (4.35),  (4.36),  (4.37), 
one  obtains  the  differential  equation  to  be  satisfied  by 

A 

[8(1  -  cos  e)  +  h  (e)o>  (e)]c  -  4sin  eg 

+  [|(1  -  3cos  e)  +  y|-h2(9H4(e)]i;  =  0  (0<|e|<n)  .  (4.39) 


The  boundary  and  transition  conditions  accompanying  (4.39)  follow  from  (4.33) 
by  recourse  to  (4.35)  to  (4.38).  In  this  manner  one  deduces 
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c(-,)43-1/4s-1/2,  c(ir)  *  ”2" 3”^ /4  ,  (4.40)1 

dO+)=dO-),  st(0+)  =  c(0-)  .  (4.41) 

We  thus  seek  a  solution  c(e)  of  the  nonlinear  boundary-value  problem  con¬ 
sisting  of  (4.39)  subject  to  (4.40) ,  (4.41 ) ,  such  that  <j>(e)  -  defined  by 
(4.37)  -  is  strictly  positive  on  [-tt,*]. 

As  far  as  the  numerical  determination  of  c(e)  is  concerned,  the  problem 
at  hand  is  awkward  because  (4.40)  prescribe  the  boundary  values  of  the  derivative 
of  the  unknown  function,  rather  than  c(±ir).  For  this  reason  we  now  convert  the 
foregoing  problem  for  g(e)  to  one  for  iji{e)  that  circumvents  this  difficulty. 
With  this  aim  in  mind  we  note  first  that  (4.37),  (4.40)  give 

*(-ir)  »  2x3_1/V1/2,  *U)  =  2x3"1/4  .  (4.42) 

We  observe  parenthetically  that  the  above  boundary  values  of  \He)  are  con¬ 
sistent  with  the  asymptotic  formulas  for  x(±ir)  arrived  at  in  (3.39),  as  is 

2 

seen,  on  setting  |a|  =  in  (3.39)  ,  by  comparing  this  estimate  with  (4.28) 
and  by  referring  to  (4.36).  Next,  we  differentiate  (4.37)  with  respect  to 
e  and  eliminate  'i  from  the  ensuing  identity  by  means  of  (4.39).  The  equation 
thus  obtained  relates  c,  i  to  i,  is  linear  in  the  former  two  variables, 
and  may  be  employed  along  with  (4.37)  to  express  c,  l  in  terms  of  ♦  ,  i  and 
trigonometric  functions  of  e.  In  this  manner  one  finds  that 


°rve  that  (4.37),  (4.38)  necessitate  c(-ir)>0  and  c(tt)<0. 

.1  that  a-j  =  0,  a^O.  See  (4.10)  and  the  remarks  following  (4.10). 


,  ■  t-  ■  in vri  i  ■  ■  i  i 
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C =  — >r - a-  ([3h^(e)iJi^  -  8cos  e  +  8]kos 

15h2(e)«r  2 

+  4[h^(e)^  -  8cos  e  +  8]isin  |-3  (-irsesit)  ,  (4.43) 

c =  — >r - tt  1C8  +  8cos  e  -3h2(e)i|j4]*sin|- 

12h2(e)4>4  2 

+  4[h2(e)\f>4  -  8C0S  6  +  8]icos  |-}  (-use**)  .  (4.44) 

Differentiating  (4.43)  and  equating  the  result  to  c  given  by  (4.44),  one  is 
led  to  the  differential  equation  for  <|»: 

[h2(o)^  +  8(l  -  cos  e)]i  -  32(1  -  cos  9)i|> 

-  4sin  Biji  -  (3  -  cos  9)iji+-^-h2(9)\jj^  =  0  ( 0  <  |  8  |  <  -rr )  .  (4.45) 

Finally,  (4.41)  together  with  (4.43),  (4.44)  furnish  the  transition  conditions 
appropriate  to  <|>: 

^(0+)  =  iJj(O-),  s^(0+)  =  i(0-)  .  (4.46) 

It  is  readily  shown  that  (4.45),  accompanied  by  conditions  (4.42),  (4.46)  and 
supplemented  by  (4.43),  In  turn  Implies  (4.39),  (4.40),  (4.41),  as  well  as 
(4.37). 

The  completion  of  the  three-term  asymptotic  representation  (4.12)  of 
y(r,e)  pertaining  to  the  special  case  b-j  =  0  has  thus  been  reduced  to  the 
task  of  finding  a  solution  «p(e )  of  the  nonlinear  second-order  differential 
equation  (4.45)  that  is  positive  for  -irsesir  and  meets  the  boundary  and 
transition  conditions  (4.42),  (4.46).  Assuming  for  the  moment  that  \p( e ) 
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has  been  so  determined,  ;(e)  is  then  supplied  by  (4.43),  z^e)  follows 
from  (4.35),  while  z2(e)  vanishes  identically  according  to  (4.34).  Further, 
the  exponents  p  and  p'  are  given  by  (4.25).  Taking  real  and  imaginary 
parts  in  (4.12)  one  therefore  has  the  following  asymptotic  results  for  the 
degenerate  case  b^  =  0  (A=0): 


y1(r,e)~t  +  cr5/,4c(e)  , 

1/2  6 

y2(r,0)~a2r/  h(0)sin  +  b2rcos  9  (-irseiir) 


>  (4.47) 


wi  th 


23'231'V'2, 


a2>0 


(4.48) 


The  first  of  (4.47)  supplants  the  weak  estimate  (4.11)  for  y-j(r,0)  furnished 
by  (4.9)  when  bj  =  0.  Note  from  (4.41)  that  c(e)  is  continuous  at  0  =  0, 
although  Its  derivative  suffers  a  finite  jump-discontinuity  there.  On  account 
of  (4.26),  (4.28),  (4.36),  and  (4.48),  d(r,e)  and  x(r,e)  at  present  conform 
to  the  estimates 


J(r,e)  ~  ^-r'1/4h(eH(e) 

(-ir  £  0  s  ir)  , 

(4.49) 

A(r,e)  ~cr1/4[h(e)i}i(0)]"1 

(-7T  S  0  £-  Tt)  , 

(4.50) 

-1/2 

whereas  (3.53)  gives  merely  the  inadequate  estimate  J(r,9)=o(r  ) 

(-itsflsn)  when  A=0.  The  result  (4.50)  reduces  to  (3.39)  for  e  =  ±rr  and 
is  consistent  with  the  a  priori  assumption  (3.40).  Since  the  right  member  in 
(4.49)  Is  the  Jacobian  determinant  of  the  right  members  in  (4.47),  the  posi¬ 
tivity  of  tjj  on  [-it ,tt]  assures  that  (4.47)  furnish  a  one-to-one  approxi¬ 
mation  to  the  local  deformation  map. 


-49- 


Our  attempts  to  deal  analytically  with  the  highly  nonlinear  two-point 
boundary- value  problem  (4.45),  (4.42),  (4.46)  ,or  ^(e)  have  remained  un¬ 
successful.  In  contrast,  no  particular  difficulties  were  encountered  In  its 
numerical  solution.  The  results  thus  obtained  leave  no  doubt  as  to  the  exist¬ 
ence  of  a  solution;  although  Its  uniqueness  cannot  be  taken  for  granted,  there 
is  no  Indication  of  non-uniqueness. 

Figure  3  and  Figure  4  display  illustrative  graphs  of  <|>(e)  and  c(e) 
for  -nsesu,  appropriate  to  two  choices  of  the  stiffness-ratio  s^i/i^* 

In  Figure  3,  s  =  1  / 2 ,  whereas  Figure  4  pertains  to  s  =  l,  that  is  to  an  all- 
around  homogeneous  Neo-Hookean  sheet.  In  the  latter  instance  both  ^(e)  and 
c(e)  are  even  functions  and  there  is  no  discontinuity  in  the  slope  of  the 
corresponding  curves  at  9  =  0. 


5.  Discussion  of  the  elastostatic  field  near  the  tip  of  the  interface-crack. 

The  special  case  s  *  1 . 

In  this  section  we  first  examine  the  structure  of  the  deformation  field 
near  the  tip  of  an  interface-crack  as  predicted  by  the  nonl  inear  theory  of  plane 
stress  for  Neo-Hookean  sheets.  Thereafter  we  shall  discuss  the  associated 
near-field  of  stress. 

For  the  purpose  at  hand  it  is  convenient  to  put  the  asymptotic  results 
(4.9)  and  (4.47)  into  a  non-dimensional  form  by  means  of  the  scaling 


2  2  2 
3  "  fc)/d2*  *”  x2^29  ^  =  r^a2 


2  2 
n-|  *  (y^  •  i)/&2*  n2  ~  ^2^d2 


>  (5.1) 


Then  (4.9),  (4.10)  yield 


'  a  j  ~  ii  i 
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n-j  (p ,0)  ~b-|pcos  6  (b^  >  0)  , 

n2(p,0) ~p^/,^h(e)s1n +  bjpcos  e  (-isesj)  , 

while.  In  the  event  that  b,  =0,  (4.47),  (4.48)  give 


y  (5.2) 


n1(p,0)~23/231/4p5/4c(e)  , 


n2(p.e)~p1,/2h(e)sin|  +  b2pcos  8  (-ttS8s.it)  . 


(5.3) 


We  observe  to  begin  with  that  the  foregoing  asymptotic  results  are  en¬ 
tirely  free  of  oscillations  of  the  kind  arising  in  the  treatment  of  the 
linearized  problem  and  discussed  in  Section  1.  Mathematically,  the  absence 
of  this  oscillatory  singular  behavior  can  be  traced  to  the  fact  that  the 
initially  undetermined  exponents  m,  m' ,  n,  n' ,  and  p,  p'  introduced  In 
(3.18),  (3.43),  and  (4.12)  were  all  found  to  be  real -valued. 

Of  particular  Interest  Is  the  shape  of  the  deformed  upper  and  lower 
crack-face  in  the  vicinity  of  the  (fixed)  end  of  the  crack.  From  (5.2)  and 


(3.34)  one  has 


n1(p.’f)','b1p,  n2(p,Tt)~p 


1/2 


n-|  (p,-if)  ~-b^p ,  ^(p.-it) — sp 


1/2 


>  (5.4) 


b.j  *  0,  s  =  /y  2  . 


j 


Elimination  of  p  between  the  first  two  and  the  second  two  of  (5.4)  leads 
to  the  subsequent  approximate  description  of  the  curves  into  which  the 
crack-faces  at  e  =  it  and  8  = -it  are  deformed: 
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ti^  *  n2 *  n2  £  ®  for  9  =  *  » 


n i  =  -  — j- n 2  »  n2  *0  for  9  =  -n 


(5.5) 


In  the  non-degenerate  case  b-|  > 0  this  description  is  adequate  and  asserts 
that  each  crack-face,  after  deformation,  is  locally  approximated  by  an  arc 
of  a  parabola,  as  indicated  in  Fig. 5.  The  two  parabolic  arcs  determined  by 
(5.5)  join  up  with  a  common  tangent  at  the  crack-tip,  this  tanqent  being 
perpendicular  to  the  crack-axis  following  the  rigid  rotation  that  led  from 
(3.51)  to  (4.9).  The  deformed  upper  and  lower  crack-faces,  when  s/1,  are 
seen  to  lie  on  two  distinct  parabolas,  both  of  which  are  concave  toward  the 
undeformed  crack. 

In  the  degenerate  case  bj  =  0,  one  draws  from  (5.3)  and  (3.34)  that 


/  i  ,3/2J/4  5/4  ,  x  /  ,  1/2 

n1(p,r)~2  3  p  cU),  n2(p.n)~P  > 

I  (p,-ir)  ~23/23ly^p5/,4c(-ir),  n2(p  ,-ir)  sp1/,?  ,  J 


(5.6) 


whence  in  first  approximation 


r»i  s  02*0  for  0  =  tt  , 

n,  -  23/231/4c(-1r)s'5/2(-n2)5/2»  n2  *  0  for  9  =  -tt  . 


>  (5.7) 


Thus  once  again  the  deformation  images  of  the  upper  and  lower  crack-faces 
have  a  common  tangent  at  the  tip  of  the  crack.  Figure  6  illustrates  the 
shape  of  the  deformed  crack-faces  when  b^  =  0,  for  the  two  stiffness-ratios 
s  =  l/2  and  s  =  l  (homogeneous  sheet).  As  is  apparent  from  the  graphs  in 
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Fig. 3  and  Fig. 4, 


c(it)  <  0,  c(-ir)  >  0  for  s  =  u-,/y2  =  l/2 


t(ir)  =  c(-tO  >  0  for  s  =  1 


r  (5.8) 


These  inequalities,  in  conjunction  with  (5.7),  imply  that  the  boundary  of 
the  sheet  near  the  tip  of  the  crack,  when  =  0,  is  deformed  into  an 
S-shaped  curve  for  ss  1/2,  whereas  for  s  =  l  the  region  occupied  by  the 
deformed  sheet  is  locally  convex.  Accordingly,  the  special  case  b^  =  0 
marks  the  transition  from  circumstances  in  which  the  deformed  sheet  is  wholly 
concave  sufficiently  close  to  the  crack-tip  to  conditions  under  which  this 
is  no  longer  true.  As  will  become  clear  later  on,  > 0  in  the  particular 
crack-problem  for  a  homogeneous  sheet  (sal)  treated  by  Wong  and  Shield  [15]. 
It  is  not  obvious,  however,  whether  or  not  there  exist  global  loading  and 
sheet  geometries  that  induce  the  degeneracy  b1  =  0. 

In  view  of  the  primary  purpose  of  this  study,  the  most  important  con¬ 
clusion  emerging  from  the  preceding  results  is  that  -  at  least  for  an  inter¬ 
face-crack  between  Neo-Hookean  sheets  -  the  finite  theory  of  plane  stress  does 
not  give  rise  to  interpenetration  of  the  deformed  crack-faces  in  the  vicinity 
of  the  crack-tips.  The  prediction  of  such  an  unacceptable  overlap  is  thus 
seen  to  stem  from  the  linearization  of  this  singular  problem.  Indeed,  some¬ 
what  surprisingly,  the  crack  is^  found  to  open  smoothly  even  if  s  1 1 . 

As  is  apparent  from  (5.2),  (5.3),  the  leading  term  in  the  corresponding 
approximations  for  n2  vanishes  when  e«0.  Consequently,  a  more  detailed 
discussion  of  the  deformations  in  the  vicinity  of  the  crack-tip  at  p  =  0 
necessitates  that  one  take  account  of  the  second  term,  whose  coefficient  b2 
in  particular  enters  into  the  lowest-order  approximation  to  the 


“’’ " v'  •"  ■  ■  *"■  ■•-  ~ '• '.-V,— ~>-  -  - 
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defonnatlon-image  of  the  bonded  interface.  Confining  this  discussion  to  the 
non-degenerate  case  b^>0,  we  gather  from  (5.2)  and  (5.1)  that 

nn  ~bi5r  r,2~H(C2)CVcf+lf  -  ^3  1/2  +  b2^  •  (5-9) 

provided  H  is  the  step-function  defined  by 

HU2)  =  1/^  if  C2  >  0,  H(C2>  =  -S//2  If  ?2<0  .  (5.10) 

Equations  (5.9)  reveal  that  the  material  coordinate  lines  =  constant 
are  —  in  first  approximation  —  carried  into  the  straight  lines  n1  =  b.|£.|, 
while  the  coordinate  lines  £2  =  constant  are  deformed  into  the  family  of 
curves  approximated  by 

n2  =  H(52)[,/(ni/bi  )2  +  cf  *  nl/,*>l3  ^  +  b2nl *  (5.11) 

Within  this  approximation  the  deformation-image  of  the  interface  S2  =  0, 

^  >  0  is  rectilinear  and  furnished  by 

^"Wbv  V0  •  (5*12) 

Figure  7  displays  qualitative  sketches,  based  on  (5.9)  to  (5.12),  of  the 
deformed  material  lines  ^  =  constant  and  ?2  =  constant,  appropriate  to 
b.j  >  0,  for  b2>0,  b2=0,  and  b2<0,  depending  on  whether  the  stiffness- 
ratio  s  <  1  or  s  =  1 . 

We  conclude  the  discussion  of  the  near-field  deformations  with  some 
observations  concerning  the  transverse  stretch-ratio  x,  which  reflects 

the  local  thinning  of  the  two  sheets.  The  asymptotic  behavior  of  x(r,e), 
as  r  +  0,  is  supplied  by  (3.62)  if  b-|  >  0  (A>0)  and  follows  from  (4.28), 
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(4.36 ) 1  when  b,  =0.  In  either  instance  A(r,e),  at  any  fixed  value  of  e 
in  [-if, if],  is  an  increasing  function  of  r  near  r=0,  so  that  the  thinning 
becomes  more  severe  as  the  crack-tip  is  approached.  Barring  the  degenerate 
case  b-|  =  0,  we  note  that  A(r,e),  for  any  fixed  (sufficiently  small)  value 
of  r,  is  smallest  at  e  =  0  and  increases  steadily  with  |e|  as  |e|-»ir. 
Consequently  the  thinning  is  most  pronounced  at  the  interface  and  less  prom¬ 
inent  at  the  crack-faces.  If  s<  1  (y1  <u,>),  the  local  thinning  is  more 

prevalent  in  the  upper  than  in  the  lower  sheet,  as  is  to  be  anticipated. 

We  turn  now  to  the  determination  of  the  stresses  near  the  crack-tip. 
Asymptotic  results  for  the  cartesian  components  of  nominal  stress  are  readily 
obtainable  from  (4.9)  with  the  aid  of  (3.13)  and  (3.40).  Thus, 


D1 


22 

U1 


=  Tj-p  1/2cos  |-  +  o(l ) 


J12 


ic  ,./i\  21- 

—  =  0(1  ),  ——  -  - 


1  -1/2 


sin- 


“2 

F[eT 


(-IT  &  0  &  1l) 


A 


Y  (5.13) 


where  p  is  the  dimensionless  radial  coordinate  introduced  in  (5.1).  The 
associated  components  of  actual  stress  are  deducible  from  (4.9)  by  means  of 
(2.44)  and  (3.40).  This  computation  yields: 


T11  .n,  -1/2,  T22  1  -lh,Q,  .  -l/2.j_ 9 

— - o(p  )»  — — "tp  h(0)-b9p  sin- 

M-J  y-j  H  L 


t12  _  t21 
yl  ’  U1 


D1  -1/2 
-g-p 


sin- 


(-t  S0  Sit) 


Y  (5.14) 


T 


See  also  Fig. 3  and  Fig. 4,  where 


*(e) 


is  plotted  for  s  =  l/2  and  s=l. 
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Evidently,  the  expansion  (4.9)  is  insufficient  to  produce  a  dominant  estimate 
for  a- 12  and  also,  when  b-j  =  0,  one  merely  infers  from  (5.13),  (5.14) 

that  o-|i  and  tend  to  zero  as  p-*0. 

Equations  (5,14)  lead  to  the  following  asymptotic  results  for  the  actual 
principal  stresses: 

~ao(p  ^^),  ~  ~  TfP  ^h(e)  -  b2P  ^sin^  (-irsesn)  .  (5.15) 

The  above  weak  estimate  for  x-|  can  easily  be  improved  if  b1  >  0.  To  this 
end  we  first  draw  from  (2.28)  that 

det  t  s  into  5  det  F  det  o  =  J  det  o  (5.16) 

and  then  Invoke  (4.10),  (3.53),  as  well  as  (5.13),  to  see  that 


1  2.  2  -1h2/„\  2  6 
TlT2~^llrlc  "  (e)cos  j  . 

On  the  other  hand,  (5.17)  and  the  second  of  (5.15)  justify 


(5.17) 


21 

U1 


~b^h(e  )cos2(e/2) 


(-TI  S  0  £  TT  ) 


(5.18) 


Since  the  deformed  crack-faces  are  free  of  tractions,  the  fact  that  t-j 
vanishes  at  e  =  ±n  is  not  accidental. 

It  Is  of  interest  to  examine  the  dependence  -  to  dominant  order  -  of 
the  primary  actual  normal  stress  X£2  upon  the  radial  distance  from  the 
crack-tip  after  deformation,  i.e.  upon  the  spatial  radial  coordinate.  We 
do  so  merely  for  an  approach  to  the  tip  along  the  Interface  e  =  0  and  for 
this  purpose  gather  from  (5.2)  that 
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P*  =  In(pfO)  |  ~  |b|p,  n=n-j+in2  • 


(5.19) 


Hence  (5.14)  now  give 

t22(p,0±)~  j-h(0±)p_1  ^hJOilo;1  ,  (5.20) 

so  that  t22(p,0)  is  both  0(p"^)  and  0(p*1)  as  p-»-0.  The  conclusion 
that  this  singularity  is  stronger  than  the  analogous  square-root  singularity 
arising  in  the  linearized  theory  is  consistent  with  earlier  findings  con¬ 
cerning  crack-tip  singularities  in  hardening  materials.^ 

We  record  next  the  asymptotic  behavior  of  the  strain-energy  density. 

By  (2.40),  (2.41), 

U  =  £(l  +  A2-3)  ,  (5.21) 

with  w=u-j  and  u=u2  in  the  upper  and  lower  half-plane,  respectively. 
Further,  from  (2.27),  (3.10),  (4.9),  and  (3.40),  one  has 

I=ya,6ya,8=  |Vy|2~iP’lhZ(e)’  X  =  0(p’1/4)  '  (5*22) 


whence 

-j-p  ^ h(e)  ,  (5.23) 

which  reflects  the  jump  discontinuity  in  U  at  9  =  0. 

The  foregoing  asymptotic  results  involve  the  three  real  amplitude 


^See  [12],  [13],  The  Neo-Hookean  material  Is  "hardening"  in  the  sense  that 
the  slope  of  the  response  curve  in  Fig. 2  is  steadily  increasing  for  \>  1. 
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parameters  a g ,  t>.| ,  b2,  the  first  of  which  -  in  view  of  (5.-)  —  governs  the 
length-scale  of  the  near-field  approximations  under  discussion.  As  has  been 
pointed  out  earlier,1  these  parameters  elude  the  local  analysis  &;>«.,  for  an 
interface-crack  between  two  semi -infinite  sheets,  are  bound  to  depend  n 
the  stiffness-ratio  s  =  vi^/u2»  the  crack-length  2i,  as  well  as  on  the  specific 
loading  at  infinity. 

The  primary  amplitude  parameter  a2  can,  in  the  usual  manner,  be  ex¬ 
pressed  in  terms  of  an  appropriate  path- independent  "^-integral"  by  recourse 
to  a  familiar  conservation  law  of  finite  elastostatics.  Let  and  C2  be 
the  two  simple  closed  curves  shown  in  Fig. 8,  which  lie  in  the  respective  half¬ 
planes  H-j  and  H2-  The  conservation  law  alluded  to  above  then  assures  that 

£(k)=J  (Unroa6n6yo  ^dA-O  (k  =  1 ,2)  ,  (5.24) 

Ck 

where  £  is  the  unit  outward  normal  vector  of  Cfc.  Next,  let  ( e )  and 
rk(e)  designate  the  rectilinear  and  circular  portions  of  Ck  labeled  in 
Fig. 8.  Noting  that  n-j  =  0  along  lk(e)  and  recalling  the  boundary  conditions 
(3.2)  together  with  the  bond  conditions  (3.4),  one  infers  that  the  contribu¬ 
tions  to  +  2^  stemming  from  L^U)  and  l2(e)  vanish,  whence 

>*l  (uV°«6V..i)d,W  •  (6'25) 

C  r(e) 

provided  C  and  r(e)  are  the  two  curves  (see  Fig. 8)  defined  by 

C«  C1  +  C2  -  L^c)  -  Lz(e)  -  r(e),  r(e)  =  I^U)  +  ^(e)  .  (5.26) 

^See  the  introductory  remarks  in  Section  4. 


SgfBB 
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Passing  to  the  limit  as  e-*-0  in  (5.25)  and  making  use  of  (5.1),  (5.2), 
(5.13),  and  (5.23),  one  readily  arrives  at 


2  8 
2  (1+s) 


2  ‘ 


(Un. 


°aeVa,l 


)di 


(5.27) 


in  which  C  is  any  simple  curve  issuing  from  an  interior  point  of  the  crack, 
terminating  at  such  a  point,  and  surrounding  the  crack-tip  situated  at 
X2  =  1. 

The  conservation  law  underlying  (5.25),  and  hence  also  this  identity, 
is  equally  valid  in  the  infinitesimal  theory  of  (generalized)  plane  stress. 
One  may  therefore  calculate  from  (5.25),  by  letting  e  tend  to  zero,  the 
value  of  ^  appropriate  to  an  available  global  solution  of  the  linearized 
interface-crack  problem  for  particular  loading  conditions.  Moreover,  this 
value  of  £  in  conjunction  with  the  first  of  (5.27)  would  yield  a  small¬ 
load  estimate  for  a2  if  one  could  take  for  granted  that  the  integrand  in 
the  second  of  (5.27)  is  approximated  uniformly  on  C  -  at  small  loads  —  by 
its  counterpart  in  the  solution  of  the  corresponding  linearized  problem. 

While  such  an  assumption  is  plausible  in  case  p-j  =  u2,  it  is  no  longer 
tenable  when  y-|  +  u2  for  curves  C  that  come  sufficiently  close  to  the 
crack-tip  since  the  elastostatic  field  predicted  by  the  linear  theory  is 
oscillatory  in  any  small  enough  neighborhood  of  the  tip,  whereas  no  such 
oscillations  occur  at  least  at  suitably  small  distances  from  the  tip  ac¬ 
cording  to  the  finite  theory.  For  this  reason,  no  matter  how  small  the 
loads,  the  linear  theory  cannot  even  supply  a  pointwise  approximation  con¬ 
sistent  with  the  finite  theory  at  all_  material  points. 

It  is  conceivable,  however,  that  despite  these  circumstances,  linear¬ 
ized  elastostatics  furnishes  a  valid  approximation  to  Q  at  small  leads. 


+.  «=-?¥S>  1^^^^te•'Hi■■=•f^r -r-  -p;  ^  -  ra^PJWn ,Tr  ~~  ”  r^ " "  ^  T  !W  T- 
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Thus,  the  global  solution  based  on  the  infinitesimal  theory  might  uniformly 
approximate  the  elastostatic  field  emerging  in  the  nonlinear  theory  on  every 
material  point  set  that  is  sufficiently  remote  from  the  tips  of  the  crack. 

The  precise  approximative  status  of  solutions  to  linearized  problems  in; - 
volving  interface-cracks  remains  an  intriguing  issue. 

Although  the  special  case  s=l  encompassed  by  the  asymptotic 

analysis  carried  out  in  this  paper,  has  been  included  in  the  preceding  dis¬ 
cussion  of  the  results  obtained,  certain  aspects  of  the  crack  problem  for  an 
all-around  homogeneous  Neo-Hookean  sheet  merit  additional  attention.  If 
s  =  1 ,  equations  (3.34),  (5.1),  (5.2)  give 


n-j  (p.Q)  ~b^p  cose  , 


n2(pie)~o^^sin|-  +  b^ pcos  e  (-isesj) 


y  (5.28) 


j 


where 

p  -  r/a2 1  n^  —  (yi_t)/®2»  ng  =  y 2/ ^2  *  ( 5 •  29 ) 

Consider  now,  in  particular,  the  Mode  I  problem  governed  by  (3.8),  so 
that  the  loading  at  infinity  is  one  corresponding  to  a  pure  homogeneous 

co 

deformation  with  the  x^-axes  as  principal  axes  and  as  principal 
stretch-ratios.  In  this  instance  one  would  anticipate  the  global  solution 
to  be  symmetric  about  the  x-j-axls.  Hence  b^  =  0  and  the  near-field  of 
deformation  predicted  by  (5,28)  bt-.~-.ss 

fr  I  n-j (p >9) ~b^p  cos  e,  n2(p»c) "-p^^sin 5-  (-isein)  .  (5.30) 

it 

This  conclusion  is  found  to  be  In  asymptotic  agreement  with  the  lowest-order 

* 

f  l 

l  i 

b-  - 
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approximation  to  the  global  solution  deduced  by  Wong  and  Shield  [15] J  pro¬ 
vided  their  results  are  adjusted  to  keep  the  crack-tip  at  x2 =  i  fixed, 
and  provided  one  takes 

^1  =  ^1  *  a2  =  •  (5.31 ) 

In  the  analogous  Mode  II  problem  (3.9)  is  in  force,  the  loading  at 
infinity  being  one  appropriate  to  a  homogeneous  deformation  of  pure  shear 

oo 

of  amount  k,  parallel  to  the  x^-axis.  The  well-known  solution  to  the 
linearized  version  of  this  problem  is  anti -symmetric  about  the  crack-axis, 
i.e.  its  displacement  field 

uo(x!  ,x2)  =  ya(xi,x2)  -  xa  (5.32) 

obeys  the  parity  relations 

u-j  (xi  ,x2)  =  -u-j  (x-| ,  -  x2) ,  u2(x-j  ,x2)  -  u2(x-|  ,  -  x2)  .  (5.33) 

Since  the  amplitude  parameters  in  (5.28)  cannot  be  chosen  so  as  to  render 
this  local  expansion  compatible  with  (5.33),  the  global  solution  to  the 
nonlinear  Mode  II  problem  at  hand  apparently  cannot  possess  the  anti-sym¬ 
metry  exhibited  by  its  counterpart  in  the  linear  theory.  This  inference  is 
also  supported  by  the  observation  that  the  nonlinear  equilibrium  equations 
(2.46),  when  cast  in  terms  of  displacements,  fall  to  be  invariant  under  the 
parity  transformation  (5.33).  Further,  whereas  the  crack-faces  fail  to 
separate  in  the  Mode  II  problem  according  to  infinitesimal  elastostatics, 
the  present  nonlinear  asymptotic  analysis  predicts  that  the  crack  does  open 

^Recall  that  the  iterative  scheme  employed  in  [15]  presupposes  large  deforma¬ 
tions  throughout  the  sheet. 
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at  least  in  a  neighborhood  of  its  ends. 

The  foregoing  conclusions  regarding  the  Mode  II  problem  for  a  Neo- 
Hookean  sheet  are  strictly  parallel  to  results  arrived  at  earlier  by 
Stephenson  [16]  in  an  asymptotic  study  pertaining  to  finite  plane  strain 
for  a  class  of  incompressible  elastic  materials.  Moreover,  Stephenson 
succeeded  in  proving  the  non-existence  of  an  anti -symmetric  global  solution 
to  this  problem  for  a  Neo-Hookean  material. 
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FIGURE  1.  INTERFACE-CRACK  BETWEEN  TWO  SEMI-INFINITE 
ELASTIC  SHEETS  AND  COORDINATES. 


FIGURE  2.  RESPONSE  OF  A  NEO-HOOKEAN 

MATERIAL  UNDER  UNI-AXIAL  STRESS 
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FIGURE  5.  DEFORMED  CRACK-FACES  FOR  b,>0 


FIGURE  6.  DEFORMED  CRACK-FACES  FOR 


MMX  - 


rp 


M 

/ 


✓  '  /  7 

/  ' 


Unclassified 


SECURITY  CLASSIFICATION  OF  This  FACE  fPTiwn  D««  Em.r.d; 


REPORT  DOCUMENTATION  PAGE 

READ  INSTRUCTIONS 

BEFORE  COMPLETING  FORM 

t.  REPORT  NUMBER  2.  OOVT  ACCESSION  NO. 

;/y  tih-Moi 

1.  recipient**  catalog  number 

4.  TITLE '  - - - —  •  - -  - - 

Large  deformations  near  the  £1ps  of  an  interface 
crack  Detween  two  Neo-Hookean  sheets •  * 

t  /  *  ~ 

».  .  TYPE  OF  REPORT  *  PERIOD  COVERED 
-/ 

«.  PERFORMING  ORO.  REPORT  NUMBER 

7.  AUTHORtAI  .  _ 

J. K. -Knowles  a-nd  El  1  /Sternberg  ;  f/iS  > 

contract  or  grant  numbs Rr«> 

N0001 4- 75-C-01 9c/ 

t.  PERFORMING  oroanization  name  ano  aoorfss 

California  Institute  of  Technology 

Pasadena,  California  91125 

10.  program  element,  project,  taik 

AREA  A  WORK  UNIT  NUMBER* 

11.  CONTROLLING  OFFICE  NAME  ANO  AOORESS 

Office  of  Naval  Research  //  . 

Washington,  D.C.  v— 

12.  REPORT  DATE 

September  *981  / 

■T*. '  NUMBER  OP-PAPEA--' 

70 

14  MONITORING  AOEnCY^AME  «  AOORE1V//  dl  Huron!  from  Con  'rolling  Oltlct) 

Q£y/2J 

1*.  SECURITY  CLASS.  fo<  Utlo  noott) 

Unclassified 

1S«.  OCCL  ASSI  FI  CATION/  DOWNGRADING 
SCHEDULE 

14.  OISTRIEUTION  STATEMENT  (ol  Ml.  Kopolt) 

Approved  for  public  release;  distribution  unlimited 

17.  Distribution  STATEMENT  (of  thm  mbsirmct  mttormd  In  Blook  20,  II  dtltmron f  from  Rm *ori) 

It.  SUPPLEMENT  ARY  NOTES 

IS.  KEY  WOROI  (Co mtlmiP  on  to  rot..  olOo  II  nocmo.orr  mr.  tdontlty  Or  Oloot  mmOot) 

Interface-crack,  elastostatics,  large  deformations 

'•  /  .  —  _  _ _  . 

SO,/  ABSTRACT  £CsnllmM  an  rorotto  tiao  II  no eaaaarr  and  Idontltr  Or  Woe*  mmOorJ 

This  paper  contains  an  asymptotic  Investigation  -  within  the  nonlinear  theory 
of  elastostatic  plane  stress  -  of  the  deformations  and  stresses  near  the  tips 
of  a  traction-free  interface-crack  between  two  dissimilar  semi-infinite  Neo- 
Hookean  sheets.  The  results  obtained  are  free  of  oscillatory  singularities 
of  the  kind  predicted  by  the  linearized  theory,  which  would  require  the  two 
deformed  faces  of  an  interface-crack  to  overlap  in  the  vicinity  of  Its  tips. 
Instead,  the  crack  Is  found  to  open  smoothly  near  its  ends,  regardless  of  the 
specific  loading  at  infinity.^ 

/  A' 


od 


roMM 
AN  7) 


1473  EDITION  OF  I  NOV  SS  p  6bSQ 
S/N  0103  LF  014-flBOI  1  ' 


.ETE 


Unclassified 


SECURITY  CL  AMI  F  1C  AT  ION  OF  THIS  RAISE  fEhoi  Oaf*  BMW 


